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ABSTRACT 

>  Resonant  laser  pumping  of  transitions  within  atoms  or  ions  represents  a 
powerful  mode  of  coupling  laser  energy  into  a  gas  or  plasma.  In  the  case  of 
laser  saturation  of  a  resonance  transition  rapid  heating  of  the  free 
electrons  by  superelastic  collision  quenching  of  the  laser  pumped  state  has 
been  shown  to  play  the  central  role  in  this  process.  We  have  developed  a 
computer  code  that  models  this  1 “laser  ionization  based  on  resonance 
saturation  (LIBORS)  and  permits  us  to  map  the  three  dimensional  nature  of 
this  complex  interaction.  We  have  been  able  to  show  by  comparison  with 
experiment  that  this  LIBORS  computer  code  can  predict  the  radial  and  axial 
electron  density  and  temperature  profiles  of  the  plasma  formed  along  the 
path  of  the  laser  pulse.  These  experimental  results  in  themselves  represent 
the  first  detailed  measurements  of  a  sodium  plasma  created  by  laser 
resonance  saturation  and  reveal  the  importance  of  laser  attenuation  in  the 
formation  of  the  plasma. 

We  have  also  shown  that  both  one  and  two  photon  resonant  laser  pumping 
of  a  cold,  unexcited  plasma  created  by  two  photon  ionization  of  strontium 
vapor  can  lead  to  rapid  excitation  of  high  lying  states  of  the  strontium 
ion.  Lastly,  we  have  demonstrated  parametric  generation  within  this  same 
strontium  plasma  when  the  laser  is  tuned  close  to  a  two  photon  transition  of 
the  ion. 


RESEARCH  OBJECTIVES 


Over  the  years  research  towards  the  development  of  short  wavelength 

(VUV  to  X-ray)  lasers  has  proceeded  down  many  paths.  Recently,  groups  at 

Livermore  and  Princeton  appear  to  have  had  some  success  using  very  powerful 
(1-31 

lasers.  '  In  each  of  these  cases  the  wavelength  of  the  pumping  laser  was 
non-specific  and  the  density  of  the  plasma  was  very  high.  Alternatively, 
the  generation  of  coherent  short  wavelength  radiation  can  be  achieved 
through  multiphoton  induced  nonlinear  processes  such  as:  two  photon 
four-wave  sum  mixing  and  third  harmonic  generation. 

The  prospects  of  achieving  high  conversion  efficiency  are  greatly 
facilitated  by  selecting  the  exciting  and  generated  fields  to  be  near 
resonant  with  allowed  optical  transitions  of  the  atom.  In  the  past,  these 
interactions  have  usually  involved  atomic  vapors, but  the  quest  for 
coherent  short  wavelength  radiation  will  invariably  require  multiphoton 
pumping  of  an  ionic  species  within  a  plasma.  However,  the  fulfilment  of 
this  promise  is  only  likely  to  be  achieved  if  a  cool  (relatively  unexcited), 
uniform  plasma  column  of  considerable  length  can  be  created.  At  modest 
densities  laser  resonance  saturation  has  been  shown  to  be  very  effective  at 
rapidly  producing  near  total  ionization  of  the  laser  pumped  species. 

The  main  objective  of  our  current  research  program  has  been  to  study 
"laser  ionization  based  on  resonance  saturation"  (LIBORS),  both 
theoretically  and  experimental ly  in  order  to  gain  a  better  understanding  of 
this  complex  laser  interaction  and  to  provide  some  insight  into  the 
capability  of  this  mode  of  ionization  to  create  long,  cool  dense  plasma 
channels.  To  accomplish  this  we  have  undertaken  the  development  of  a 
computer  code  that  attempts  to  model  LIBORS  in  both  space  and  time.  We  have 
also  performed  a  set  of  experiments  designed  to  evaluate  the  properties  of 
the  plasma  channel  created  by  LIBORS  and  to  test  the  predictions  of  our 
computational  model  in  this  regard. 

Laser  saturation  of  an  ion  resonance  line  within  a  plasma  offers  the 
prospect  of  very  rapid  heating  of  the  free  electrons  through  superelastic 
collision  quenching  of  the  laser  pumped  resonance  state  population. This 
mode  of  coupling  laser  energy  into  a  plasma  could  be  very  competitive  with 
inverse  bremss4  ahlung^17^  and  could  lead  to  a  form  of  shock  heating  of  the 
free  electrons.  If  the  initial  plasma  was  cold  and  relatively  unexcited,  it 


might  be  possible  with  this  technique  to  create  a  population  inversion  that 
could  lead  to  efficient  short  wavelength  laser  action.  Consequently, 
another  of  our  objectives  has  been  to  undertake  a  preliminary  investigation 
of  laser  pumping  of  an  ion  resonance  line  within  a  plasma. 

Near  resonance  multiphoton  excitation  can  also  lead  to  certain  enhanced 
nonlinear  interactions  and  we  have  undertaken  a  preliminary  investigation  of 
one  such  process. 


STATUS  REPORT 

Laser  Ionization  Based  on  Resonance  Saturation  (LIBORS)  Program 

Laser  saturation  of  a  resonance  transition  within  an  atom  or  an  ion 
represents  a  new  class  of  laser-plasma  interactions  In  the  case  of  a 
neutral  atomic  vapor,  laser  resonance  saturation  gives  rise  to  many 
processes  that  quickly  lead  to  ionization. *la)  Once  an  appreciably  degree 
of  ionization  has  been  achieved  the  primary  mechanism  for  converting  laser 
energy  into  plasma  energy  is  superelastic  electron  collision  quenching  of 
the  laser  pumped  resonance  state  population.  This  is  an  extremely  efficient 
method  of  laser  energy  transfer  which,  in  effect,  rap.dly  cycles  the  bound 
electrons  between  the  ground  and  resonance  states,  converting  the  energy  of 
one  photon  into  free  electron  kinetic  energy  during  each  cycle  per 
atom.*9”11) 

Laser  resonance  saturation  can  be  viewed  as  complementary  to  inverse 
bremsstrahl ung ,  for  although  it  can  lead  to  energy  deposition  that  is  many 
orders  of  magnitude  greater  than  attainable  through  inverse  bremsstrahl ung 
for  modest  values  of  laser  irradiance  —  the  free  electron  temperature 
achieved  is  limited  by  the  characteristics  of  the  atom  or  ion  being 


i rradiated. 


Laser  saturation  of  an  ion  resonance  line  within  a  plasma 


saturation  can  lead  to  extremely  rapid  electron  heating,  in  effect,  creatiny 
a  well  defined  electron  temperature  jump*1**)  and  a  commensurate  sharp  burst 
in  the  emission  of  short  wavelength  radiation  from  a  suitable  plasma. 

The  objective  of  our  research  program  is  to  study  this  new  form  of 
laser-plasma  interaction  and  explore  some  of  its  potential  applications.  In 
any  practical  application,  the  laser  pulse  has  to  propagate  through  a  finite 
depth  of  the  medium  being  irradiated  and  as  a  result  of  the  strong  nature  of 
the  resonance  interaction,  the  laser  pulse  will  undergo  an  appreciable 


degree  of  attenuation  and  temporal  and  spatial  distortion.  Consequently, 
any  theoretical  modeling  of  this  interaction  should  take  account  of  this 

fact. 

To  date,  we  are  the  first  to  develop  a  computational  analysis  that 
attempts  to  model  the  three  dimensional  nature  of  this  interaction  by  taking 
account  of  the  attenuation  and  distortion  suffered  by  the  laser  pulse  as  it 
burns  its  way  through  a  realistic  atomic  vapor  distribution.  In  order  to 
assess  the  reliability  of  our  theoretical  work  we  are  comparing  our  computer 
simulations  with  the  results  attained  from  experiments  involving  a  new 
facil ilty. 

LIBORS  Experimental  Facility 

This  new  facility  basically  comprises:  a  Nd-YAG  laser  pumped  dye 
laser,  a  specially  designed  sodium  heat  sandwich  oven  which  provides  360° 
optical  access  and  a  photodetection  system  tht  includes  an  RCA  C31034 
photomultiplier  mounted  on  the  exit  slit  of  a  SPEX  monochromator.  A 
photograph  of  the  facility  is  displayed  as  figure  1,  and  a  schematic 
overview  is  presented  in  figure  2.  Detailed  descriptions  of  this  facility 
are  provided  in  Appendices  A,  B  and  C  —  which  constitute  three  papers 
describing  our  experimental  work.  Appendix  A  is  a  paper  that  has  recently 
been  published,  while  Appendices  B  and  C  represent  two  new  papers  that  will 
shortly  be  submitted  for  publication. 

The  sodium  atom  density  distribution  across  this  novel  oven  was 
evaluated  by  measuring  the  absorption  of  a  narrow,  well  collimated  beam  of 
continuum  radiation  from  a  xenon  lamp.  These  measurements  were  undertaken 
at  a  number  of  locations  across  the  heat  sandwich  oven  and  are  made  possible 
by  the  design  of  this  oven.  The  width  of  the  spectral  hole  created  by  the 
sodium  resonance  line  was  evaluated  for  the  beam  passing  through  the  high 
density  vapor  core,  while  the  depth  of  the  589  nm  spectral  hole  was 
determined  in  the  low  density  rim  of  the  vapor  disc.  The  experimental 
technique  is  described  in  detail  in  Appendix  A. 

Our  principal  plasma  measurements  have  involved  evaluating  the  radial 
profiles  of  the  free  electron  density  and  temperature  along  the  path  of  the 
laser  pulse.  We  have  also  looked  with  some  depth  into  the  accuracy  of  such 
measurements  and  the  care  needed  in  their  interpretation.  This  work  is 
presented  briefly  in  the  body  of  this  report  as  the  full  details  are 


included  in  Appendices  B  and  C  in  the  form  of  papers  that  will  shortly  be 
submitted  for  publication. 


Electron  Density  Radial  Profile  Measurements 

We  have  undertaken  the  first  spatially  resolved  measurement  of  the  free 
electron  density  in  a  sodium  plasma  produced  by  laser  resonance  saturation. 
This  evaluation  involved  fitting  the  radially  inverted  experimental  42D-32P 
multiplet  spectrum  with  convoluted  Stark  (electron  impact)  and  instrument 
broadened  profiles.  A  pair  of  representative  radial  profiles  of  the  free 
electron  density,  evaluated  at  two  equidensity  locations  along  the  path  of 
the  laser  pulse  (z  «  -2  and  2  cm)  are  displayed  as  figure  3.  These  results 
clearly  reveal  that  the  plasma  narrows  and  the  peak  degree  of  ionization 
diminishes  with  axial  location.  We  attribute  this  to  absorption  of  the 
laser  pulse  as  it  propagates  through  the  sodium  vapor.  A  detailed 

(19) 

description  of  this  work  is  presented  as  Appendix  B  and  Cappelli. 

Electron  Temperature  Measurements 

The  free  electron  temperature,  Te,  of  a  sodium  plasma  created  by 

laser  resonance  saturation  has  been  evaluated  by  two  independent  methods. 

The  first  involves  the  measurement  of  the  spectrally  integrated  emission  of 

several  sodium  multiplets  and  the  deduction  of  the  Te  from  the  slope  of 

the  least  square  linear  fit  to  this  data  when  plotted  logarithmically  (this 

is  termed  a  Boltzmann  plot).  A  representative  example  of  such  a  measurement 

for  our  plasma  is  presented  as  figure  4.  Although  this  technique  has  been 

(20  21) 

used  by  a  number  of  researchers/  *  '  we  have  shown  (see  Appendix  C)  that 

this  kind  of  measurement  can  be  subject  to  appreciable  errors  if  the  plasma 
is  not  uniform  and  optically  thin. 

In  particular,  we  have  undertaken  a  numerical  simulation  of  the 

emission  observed  from  a  cylindrical  plasma  and  have  shown  that  if  the  core 

is  highly  ionized,  then  the  bulk  of  the  emission  originates  from  a  cooler 

outer  angular  zone  where  the  upper  state  population  density  is  a  maximum 

and,  under  these  circumstances,  the  measured  temperature  fails  to  reflect 

the  value  in  the  plasma  core.  This  could  explain  some  of  the  low 

(?o  211 

temperatures  observed.'  ’  '  In  addition  self  absorption  has  been  found  to 

lead  to  a  large  scatter  of  the  multiplet  emission  data  on  a  Boltzmann  plot. 


and  even  more  deviation  from  linearity  can  arise  in  high  density  plasmas 
(where  Stark  broadening  can  be  significant)  if  the  photodetection  bandwidth 
is  inadequate. 

The  second  technique  is  indirect,  nevertheless  it  provides  us  with  the 
radial  variation  of  Te.  In  this  approach  the  free  electron  density  radial 
profile,  discussed  above  and  in  Appendix  B,  is  used  in  conjunction  with  the 
sodium  atom  density  and  the  Saha-equation ,  assuming  the  establishment  of 
Local  Thermodynamic  Equilibrium  (LTE),  to  calculate  the  free  electron 
temperature  radial  profile  T  .  Two  such  T  -radial  profiles  corresponding  to 
the  two  Ne-radial  profiles  displayed  in  figure  3,  are  presented  in  figure 
5  and  reveal  that  the  temperature  tends  to  decrease  monitonical ly  with 

increasing  radius  and  also  decreases  along  the  path  of  the  37  mJ  laser  pulse 
(from  z  =  -2  to  2  an).  The  peak  sodium  density  in  this  run  was  also 

1.5*1016  cm-3  and  the  value  at  the  axial  locations  of  interest  was  about 

1016  cm*3.  Furthermore,  since  these  measurements  were  undertaken  about  65 
ns  after  the  start  of  the  laser  pulse  (duration  40  ns)  the  plasma  is 

expected  to  be  in  LTE. 

The  observed  decline  of  the  free  electron  temperature  along  the  path  of 

the  laser  pulse  represents  an  important  observation  for  it  also  reflects 

the  progressive  attenuation  of  the  laser  pulse  with  increasing  penetration 

of  the  sodium  vapor.  Moreover,  the  low  value  of  the  electron  temperature 

corresponding  to  propagation  distances  of  several  centimeters  at  a  mean 

sodium  atom  density  of  around  1016  cm*3  is  similar  to  that  observed  by 

(20  211 

others  under  similar  experimental  conditions.  *  ' 


Multishot  Average  Multiplet  Spectrum 


In  our  experimental  measurements,  as  with  many  others,  the  spectral 
distribution  of  the  emission  is  determined  from  a  series  of  measurements  at 
different  wavelengths.  If  the  shot-to-shot  reproducibility  is  poor,  or  the 
signal  is  weak,  the  emission  recorded  in  each  spectral  interval  is  an 
average  of  several  shots.  When  the  results  of  such  experiments  are  to  be 
compared  with  those  predicted  by  a  computational  model  the  question  of 
whether  the  laser  pulse  can  be  represented  by  a  Gaussian  radial  profile 
arises.  This  is  particularly  true  if  the  interaction  is  highly  nonlinear, 


as  in  the  case  of  LIBURS,  ana  if  the  experimental  laser  pu'.  st  :.as  a  noisy 

radial  profile. 

We  have  undertaken  a  computational  simulation  of  this  problem  and  have 
shown,  in  Appendix  D,  that  computer  modeling  of  LIBORS  based  on  a  Gaussian 
approximation  for  the  laser  pulse  radial  profile  is  acceptable  provided  the 
experimental  results  are  averaged  over  at  least  16  shots. 

LIBORS  Computer  Code 

We  have  developed  the  first  computer  code  for  modeling  the  formation  of 
a  cylindrical  sodium  plasma,  through  laser  resonance  saturation,  that  takes 
account  of  the  nonuniform  sodium  atom  distribution  along  the  path  of  the 
laser  pulse  and  the  temporal  and  spatial  distortion  of  the  laser  pulse 
arising  from  absorption  as  it  propagates  through  the  sodium  vapor. 

This  LIBORS  computer  code  represents  a  major  development  of  our 

(22) 

research  group  and  is  described  in  detail  in  Appendix  E  and  Wong  and 

(o-i  \ 

Kissack.'  '  The  predictions  of  this  code  have  provided  us  with  the  first 
insight  into  the  three  dimensional  nature  of  LIBORS  and  have  been  tested  by 
comparison  with  the  results  of  experiments  using  the  facility  described 
above  and  in  Appendices  A,  B,  C  and  E.  One  of  the  most  significant  results 
of  this  work  has  been  the  recognition  that  considerable  attenuation  of  the 
laser  pulse  can  arise  within  a  few  centimeters  of  vapor  at  a  density  of 
about  1016  cm-3  and  that  this  can  lead  to  the  formation  of  a  low  temperature 
low  density  plasma  at  these  locations. 

An  idea  of  the  predictive  capability  of  this  code  is  provided  in 
figures  6  and  7  which  reveal  two  views  of  the  plasma  channel  expected  to  be 
formed  in  one  of  our  experiments.  In  figure  6,  the  axial  variation  of  the 
electron  density  for  several  radial  positions  is  displayed,  while  in  figure 
7  the  radial  variation  of  the  electron  density  is  presented  for  several 
axial  locations.  Also  shown  in  figure  6  is  the  experimentally  based  sodium 
atom  density  distribution  along  the  path  of  the  laser  pulse  used  in  the 
computer  code. 


Computational  and  Experimental  Comparison 

In  order  to  test  the  modeling  capability  of  our  LIBORS  computer  code  we 
have  compared  its  predictions  regarding  the  radial  and  axial  distributions 


of  the  free  electron  density  and  temperature  with  the  results  obtained  from 
an  experiment.  In  general  the  agreement  is  very  good  when  the  uncertainty 
in  the  experimental  measurements  is  taken  into  account  and  it  is  borne  in 
mind  that  many  of  the  cross-sections  used  in  the  code  are  not  accurately 
known.  The  20%  uncertainty  in  the  sodium  atom  density  measurements  are 
particularly  troublesome  because  this  distribution  is  used  in  the  computer 
code  and  as  a  consequence  leads  to  two  kinds  of  theoretical  uncertainty:  a 
direct  uncertainty  in  the  local  value  of  the  free  electron  density  and  an 
uncertainty  that  results  from  the  cumulative  uncertainty  in  the  laser  energy 
fluence  that  penetrates  to  that  location. 

With  this  in  mind  we  present  in  figure  8  a  comparison  of  the 
experimentally  measured  (see  Appendix  B)  and  computer  predicted  electron 
density  radial  profile  at  Z  =  -2  cm  and  t  =  65  ns  for  the  sodium  atom 
distribution  displayed  in  figure  6.  The  two  computed  profiles  are  seen  to 
flank  the  experimental  values  and  arise  from  using  two  values  for  the 
cross-section  of  the  dominant  seed  ionization  (laser  assisted  associative 
ionization).  The  computer  code  based  on  the  cross  section  derived  by 
Roussel  et  al  ,  1980,  a A  =  a^(  F)  [5x10" 15  cm2  for  most  of  our  conditions] 
fits  closest  to  the  experimental  data  on  axis,  while  for  large  radial 
positions  this  code  overestimates  the  electron  density.  Reference  to  figure 
9  shows  good  agreement  between  the  predictions  of  this  code  and  the 
experimentally  determined  axial  profiles  of  the  free  electron  density  and 
temperature.  We  believe  that  the  divergence  observed  for  large  Z  could 
result  from  error  in  the  sodium  atom  density  measurements  while  the  radial 
difference  could  be  attributed  to  an  overestimate  in  the  radius  of  the  laser 
beam.  The  experimental  electron  temperature  data  were  derived  from  the 
electron  density  data  assuming  LTE.  A  more  detailed  discussion  of  this  work 
is  provided  in  Appendix  E,  which  will  serve  as  the  basis  of  a  paper  to  be 
submitted  for  publication. 


Development  of  New  LIBORS  Computer  Code  Including  Thermal  Conduction  and 


Using  an  Adaptive  Grid 


The  previous  section  clearly  indicates  that  our  current  LIBORS  computer 
code  is  capable  of  predicting  with  reasonable  accuracy  the  long  term  state 


of  a  plasma  created  by  laser  resonance  saturation.  However,  for  certain 


kinds  of  application  the  state  of  the  plasma  during  the  formation  stage 
could  be  of  interest.  Indeed,  some  applications  may  even  rely  on  the  kind 
of  temperature  shocks  predicted  to  occur  during  the  formation  phase  of  the 
plasma.  For  example  we  present  in  figure  10(a)  and  (b)  the  radial  variation 
of  the  free  electron  density  and  temperature  at  t  =  58  ns  and  Z  *  0  where 
N0  =  1.7xl016  cm-3.  The  incident  laser  pulse  was  assumed  to  have  an 
energy  of  37  mJ  and  its  radial  distribution  of  energy  fluence  is  also 
depicted  in  figure  10(a).  In  figure  10(b)  we  display  the  corresponding 
radial  variation  of  the  thermal  conduction  and  net  time  rate  of  change  of 
the  free  electron  thermal  energy  density.  This  figure  clearly  indicates 
that  thermal  conduction  can  become  very  large  and  even  dominate  in  the 
energy  equation  where  steep  gradients  in  Te  exist. 

We  surmize  that  the  omission  of  thermal  conduction  and  electron 

diffusion  may  not  be  justified  during  the  formation  phase  and  we  are  in  the 

process  of  including  these  effects  into  the  energy  and  ionization  equations. 

These  major  advances  in  our  theoretical  treatment  have  forced  us  to 

reformulate  our  LIBORS  computer  code  so  that  integration  in  time  is  now 

undertaken  at  a  set  of  so  called  "nodal  positions"  in  space.  To  accomplish 

this  with  a  minimum  use  of  computer  time,  we  are  employing  an  "adaptive 

grid"  to  constantly  move  a  limited  number  of  nodal  positions  into  the 

regions  of  strong  gradients.  In  this  way  the  code  should  be  able  to 

accommodate  the  formation  and  propagation  of  the  ionization  front  and  the 
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associated  steep  gradient  in  the  electron  temperature,  Kissack.  ' 


MULTIPHOTON  INTERACTIONS  IN  STRONTIUM 
Strontium  Experimental  Facility 

Our  experimental  facility  for  studying  the  formation  of  a  strontium 
plasma  and  the  subsequent  one  and  two  photon  excitation  of  the  SrII  ions  is 
schematically  presented  as  figure  11.  The  strontium  vapor  is  contained  in  a 
specially  designed  crossed  heat  pipe  heated  to  around  800°C.  Our  Nd-YAG 
laser  (JK  HY750)  is  used  with  a  DCDA  doubler  and  KD*P  mixer  crystal  to 
generate  third  harmonic  radiation  at  355  nm  for  pumping  Stillbene  420  within 
our  dye  laser  (Quanta-Ray  PDL-1). 

The  performance  of  this  dye  declines  fairly  rapidly,  the  laser  output 
dropping  by  about  50%  within  20  minutes.  The  subsequent  decline  is  then 


much  slower.  For  most  of  our  initial  work  we  have  used  an  output  pulse 
energy  of  about  2-5  mJ.  A  lens  was  used  just  ahead  of  the  crossed  heat  pipe 
to  focus  the  laser  pulse  and  produce  maximum  two  photon  ionization  in  the 
centre  of  the  oven  in  order  that  the  side  light  fluorescence  could  be 
observed. 

A  SPEX  1700  II  monochromator  was  employed  to  spectrally  select  the - 
radiation  detected  by  our  RCA-4526  photomultiplier.  Radiation  emitted 
col  inear  and  parallel  with  the  laser  pulse  was  sampled  by  a  beam  splitter 
and  two  mirrors.  The  mirror  positioned  in  front  of  the  monochromator  is 
rotatable  through  90  degrees  so  that  it  can  also  be  used  in  conjunction  with 
another  beam  splitter  and  mirror  to  sample  any  radiation  emitted  in  the 
opposite  direction  to  the  laser  pulse.  This  mirror  can  also  be  removed  and 
a  lens  introduced  to  focus  some  of  the  strontium  plasma  side  light  emission 
(i.e.,  that  emitted  perpendicular  to  the  laser  beam)  onto  the  entrance  slit 
of  the  monochromator.  Two  photodiodes  were  also  used  in  the  experiments: 
one  to  monitor  the  output  of  the  dye  laser,  the  other  was  used  to  provide  a 
trigger  pulse  to  the  signal  averager  (EG&G  4420)  used  for  data  acquisition 
of  the  photomultiplier  signal. 


Superelastic  Heating  of  a  Strontium  Plasma  by  Resonant  One  and  Two  Photon 


A  number  of  years  ago  we  had  theoretically  predicted  that  laser 
saturation  of  an  ion  resonance  line,  within  a  dense  plasma,  could  lead  to 
extremely  rapid  rates  of  electron  heating  through  superelastic  collision 
quenching  of  the  laser  pumped  state. We  had  also  shown  that,  under 
certain  circumstances,  this  superelastic  laser  energy  conversion  process 
could  couple  laser  energy  into  the  plasma  several  orders  of  magnitude  faster 
than  through  inverse  bremsstrahl ung  .  ^ 1 7 )  Furthermore,  we  suggested  that 
multiphoton  laser  pumping  might  be  used  to  expand  the  capabilities  of  this 
interaction.^17) 

We  have  recently  undertaken  the  first  test  of  superelastic  electron 
heating  of  a  strontium  plasma  created  by  two  photon  ionization  of  strontium 
vapor.  The  details  of  this  work  are  provided  in  Appendix  F,  which 
constitutes  a  paper  that  will  shortly  be  submitted  for  publication  and 
Menard.  ^*) 


In  essence  strontium  vapor  confined  within  our  crossed  heat  pipe  oven 
is  two  photon  ionized  with  a  pulse  of  laser  radiation  of  wavelenyth  (^) 
less  than  431  nm.  When  the  wavelenyth  of  the  laser  was  tuned  to  the 
52P1/2-52S1/2  SrII  resonance  transition  at  421.67  nm,  strong  emission  was 
observed  from  the  high  lying  transition,  52F5/2, 7/2-s2l)5/2»  at  nm. 

We  believe  that  this  excitation  represents  clear  evidence  of  free  electron 
heating  resulting  from  superelastic  electron  collision  quenching  of  the 

laser  pumped  52P1/2-52S1/2  transition.  The  52S1/2  ion  ground  state 
population  being  created  by  two  photon  ionization  with  the  same  laser  pulse. 
We  were  able  to  eliminate  the  possibility  that  this  565.2  nm  emission  arose 
through  three  photon  excitation  on  the  grounds  that  the  intensity  of  the 
emission  failed  to  peak  as  the  laser  was  tuned  through  the  wavelength  for 

exact  three  photon  resonance  at  422.14  nm,  while  it  did  peak  when  the  laser 

was  tuned  through  the  421.67  nm  wavelength  corresponding  to  the  ion 
resonance  wavelength,  see  figure  12(a).  Although  the  emission  increases 
sharply  as  the  laser  wavelength  approaches  421.67  nm  from  the  blue,  the 
emission  falls  much  more  slowly  on  the  red  side  of  421.67  nm.  We  attribute 
this  asymmetry  to  some  form  of  sel f- focussing  effect.  Both  the  intensity  of 
this  emission  and  its  spectral  asymmetry  are  seen  to  decline  with  decreasing 
strontium  density,  figure  12.  The  most  convincing  evidence  that  the 
excitation  mechanism  is  collisional  in  origin  rather  than  laser  induced 
stems  from  the  long  lived  nature  of  the  emission. 

The  emission  at  565.2  nm  is  also  found  to  peak  sharply  when  the  laser 
is  tuned  through  the  two  photon  52S1/2-62S1/2  transition  at  418.97  nm,  which 
leads  us  to  suspect  that  again  superelastic  collision  quenching  of  the  laser 
pumped  62S1/2  level  population  could  be  heating  the  free  electrons  and 
enabling  them  to  excite  the  high  lying  52F  levels.  In  this  instance  the 
565.2  nm  emission  arises  only  when  the  laser  wavelength  exactly  matches  the 
value  needed  to  excite  the  two  photon  transition  because  only  then  is  the 
62S1/2  state  actually  populated.  Verification  of  the  collisional  origin  of 
the  interaction  is  again  provided  by  the  long  lived  nature  of  the  565.2 


emission. 


We  conclude  that  this  work  constitutes  the  first  proof  of  principle 
that  one  or  two  photon  laser  pumping  of  an  ion  can  lead  to  superelastic 
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heating  of  the  free  electrons  and  concomitant  rapid  ex^itction  of  high  lying 
ion  states. 

Parametric  Generation  Within  a  Strontium  Plasma  Created  by  Two-Photon 
Ioni zation 

The  study  of  nonlinear  optical  effects  in  plasmas  has  in  the  past - 
primarily  been  associated  with  the  nonlinear  polarization  of  the  free 
electrons.  However,  when  the  exciting  laser  is  tuned  to  be  in  close 
resonance  with  a  two  photon  transition  within  the  ions  of  the  plasma  the 
nonlinear  susceptibility  of  the  ions  can  dominate.  As  part  of  our  study  of 
resonant  multiphoton  interactions  in  strontium  we  have  found  that  an  intense 
beam  of  radiation  almost  col  inear  with  the  pump  laser  beam  can  be  generated 
parametrically  within  a  strontium  plasma  that  is  created  by  two  photon 
ionization.  What  is  particularly  interesting  is  that  this  four-wave 
difference  mixing  is  achieved  with  the  same  laser  that  is  responsible  for 
the  ionization  of  the  strontium. 

This  complex  sequence  of  interactions  can  be  understood  by  reference  to 
figure  13  which  shows  how  laser  radiation  at  418.97  nm  can  two  photon  ionize 
Sri  and  excite  the  four-wave  difference  mixing  process  leading  to  beams  of 
radiation  at  430.67,  416.30,  407.89  and  421.67  nm  that  are  close  to  col  inear 
with  the  laser  pump  beam.  The  details  of  this  experiment  and  a  formulation 
of  its  theoretical  basis  is  provided  as  Appendix  G,  which  constitutes 
another  paper  that  will  shortly  be  submitted  for  publication,  and  in  the 
M.A.Sc.  thesis  of  K.  Menard. 
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NEW  DISCOVERIES  STEMMING  FROM  RESEARCH 

Saturation  of  an  atomic  transition  by  the  intense  radiation  field  of  a 
suitably  tuned  laser  represents  an  important  kind  of  interaction  with  a  wide 
range  of  potential  applications.  The  consequences  of  laser  resonance 
saturation  and  the  applications  stemming  from  this  interaction  depend  to  a 
very  large  extent  upon  the  period  of  saturation.  If  the  resonance  to  ground 
level  populations  are  only  momentarily  locked  in  the  ratio  of  respective 
degeneracies  (saturation)  the  principal  effect  is  a  burst  of  intensified 

t  ]  ?) 

spontaneous  emission  that  can  be  used  to  diagnose  the  excited  medium.'  *  ; 
On  the  other  hand  an  extended  period  of  saturation  (lasting  for  much  longer 
than  the  resonance  state  lifetime)  can  lead  to  extensive  perturbation  of  the 
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mtu  i  uni .  ’  Indeed,  if  the  free  electron  superelastic  collision  time  is 

short  compared  to  the  duration  of  .the  laser  pulse  near  total  ionization  of 

(4-19 1 

the  laser  pumped  species  is  rapidly  achieved.  '  In  the  case  of  an  ionic 
species  extremely  rapid  changes  of  electron  temperature  can  be  produced  by 
this  means. 

The  author  was  the  first  to  recognize  the  importance  of  laser 
saturation  and  many  of  its  possible  areas  of  appl ication.^"^  Subsequent 
work  in  the  author's  laboratory  and  elsewhere  have  proven  that  laser 
resonance  saturation  does  represent  a  significant  form  of  interaction 
between  laser  radiation  and  atomic  vapors  or  plasmas. 

Momentary  laser  saturation  represents  a  powerful  diagnostic  technique 

that  is  finding  application  in  many  areas  ranging  from  fusion  reactor 

(211  (22  231 

studies'  '  to  combustion  measurements.'  *  '  The  combination  of  laser 

ablation  and  laser  saturation  spectroscopy  represents  a  new  approach  at 

evaluating  fundamental  atomic  quantities  such  as:  radiative  lifetimes, 

branching  ratios,  transition  probabilities  and  selected  collision 

cross-sections.  A  preliminar  paper  on  this  subject  was  published  by  us  in 
(241 

Physical  Review.'  '  This  technique,  as  well  as  being  convenient  and 
accurate,  is  particuarly  well  suited  for  measurements  on  short  lived,  highly 
ionized  species  created  by  laser  ablation.  Furthermore,  it  is  versatile  and 
can  use  multiphoton  or  stepwise  excitation  as  the  means  of  generating  the 
bursts  of  intensified  emission. 

As  a  spin-off  of  this  work  we  have  also  shown  that  this  concept  can 
also  form  the  basis  of  a  new  form  of  trace  element  laser  microprobe  called  a 
TABLASER.(25) 

More  recently,  we  have  demonstrated  in  a  proof  of  principle  experiment 
that  laser  saturation  spectroscopy  can  be  used  to  directly  measure,  with 
both  spatial  and  temporal  resolution,  the  ion  to  neutral  atom  density  ratio 
in  a  rapidly  expanding  ablation  plasma. 

The  ionization  capabilities  of  extended  laser  resonance  saturation  have 
now  been  indisputably  demonstrated  with  a  variety  of  experiments. ye 
have  developed  a  model  of  this  "laser  ionization  based  on  resonance 
saturation"  (LIBORS)  process  and  have  shown  that  the  interaction  can  be 
thought  to  proceed  in  four  stages. Confirmation  of  certain  aspects 
of  our  theory  has  been  provided  by  several  research  groups. Recently, 
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there  has  been  some  success  in  usiny  laser  resonance  saturation  for  the 

(27) 

purpose  of  creating  infrared  laser  action.  ' 

Recently  we  have  developed  a  new  form  of  the  LIBORS  computer  code  that 
enables  us  to  gain  some  insight  into  the  three  dimensional  nature  of  this 
interaction.  The  predictions  of  this  code  have  been  found  to  be  in 
reasonable  agreement  with  experimental  measurements  of  the  radial  and  axial 
electron  density  and  temperature  of  the  cylindrical  plasma  created  by  laser 
resonance  saturation.  LIBORS  is  currently  viewed  as  one  of  the  least  risky 
methods  of  producing  the  cold  lithium  plasma  used  for  ion  beam  generation  in 
the  Sandia  National  Laboratory's  PBFA  II. 

Recently,  we  have  demonstrated  that  one  and  two  photon  resonant  laser 
pumping  of  cold,  unexcited  ions  formed  by  two  photon  ionization  can  lead  to 
rapid  excitation  of  high  lying  states  of  these  ions  (Appendix  F).  We 
believe  that  this  approach  could  be  used  to  produce  intense  sources  of  short 
wavelength  radiation  and  possibly  an  efficient  short  wavelength  laser.  We 
have  also  shown  that  two  photon  resonant  pumping  of  the  62Si/2-52S1/2 
forbidden  transition  of  SrII  in  a  cold  strontium  plasma  can  lead  to 
parametric  generation  of  intense  radiation  that  is  almost  col  inear  with  the 
pump  laser  beam. 
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Figure  3.  Free  Electron  Density  Radial  Profiles  Attained  at  Z  =  -2  and  2  cm 
for  N0  =  10locm-3  and  t  =  65ns,  E^=25mJ. 
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Logarithmic  variation  of  the  observed  emission  intensity  of  several 
sodium  multiplets  with  the  energy  of  excitation.  Also  shown  is 
the  least  square  linear  fit  to  the  calculated  data  -  the  slope  of 
this  line  gives  an  electron  temperature  of  5496  K. 
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Fig.  5  Representative  Free  Electron  Temperature  Radial  Prof iles Deduced 
From  the  Measured  Electron  Density  Radial  Prof iles  Assuming  LTE. 
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Fig.  8  Free  Electron  Density  Radial  Profile  at  Z  =  -2cm 

(F)al^^e^nSy<;A  --XPeriment  (#) 


a 

8 


'AS 


Glass 


Mirror 


N0  =  2.8  x  lO^cm'* 


TUNING  WAVELENGTH 


Fig.  12 


Variation  of  Sr  II  ( 52F5/2 ,7/^“5^Dc/2)  565.20nm  Emission  from 
a  Strontium  Plasma  as  a  Function  of  the  Exciting  Laser  Wave¬ 
length  for  Three  Values  of  the  Initial  Strontium  Atom  Density. 
Ion  Resonance  Line  Wavelength  at  421.67nm. 
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Fig.  13  Partial  Grotrlan  energy  level  diagram  of  Sri  and  II, 
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Sodium  atom  distribution  within  a  heat  sandwich  oven 
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In  a  heat  sandwich  oven  the  metal  vapor  is  confined  to  be  disk  shaped  and  is  optically  accessible 
through  360°  in  the  plane  of  the  disk.  We  have  used  this  feature  to  optically  measure  the  radial 
atom  density  distribution  of  sodium  vapor  within  this  type  of  oven,  under  a  range  of  conditions.  In 
particular,  we  have  observed  the  formation  of  a  donut-shaped  atom  density  distribution  when  the 
heat  sandwich  oven  is  operated  at  high  temperatures  with  an  under  pressure  of  the  argon  buffer 
gas. 


INTRODUCTION 

Alkali  metal  vapors  are  highly  reactive  and,  in  the  absense  of 
adequate  precautions,  can  quickly  spoil  'he  windows  of  any 
containment  vessel  built  for  optically  studying  them.  One 
method  of  minimizing  this  degradation  is  to  interpose  a 
cooled  inert  buffer  gas  between  the  hot  vapor  and  the  win¬ 
dows.  Vidal  and  Cooper1  first  demonstrated  that  a  combina¬ 
tion  of  the  heat  pipe  concept  with  a  buffer  gas  not  only  pro¬ 
tected  the  windows  but  led  to  the  formation  of  a  well-defined 
column  of  the  vapor. 

The  limited  optical  access  of  this  system  was  overcome 
by  Boyd  et  at 2  who  proposed  a  sandwich-like  oven  in  which 
the  vapor  was  trapped  between  two  flat  disk-shaped  wire 
mesh  wicks  that  were  heated  in  their  center  and  cooled 
around  their  periphery.  This  design  leads  to  the  creation  of  a 
disk-shaped  region  within  the  oven  and  to  360°  optical  access 
in  the  plane  of  the  disk. 

I.  HEAT  SANDWICH  OVEN 

We  have  designed  and  built  our  own  heat  sandwich  oven 
for  undertaking  experiments  involving  laser  ionization  of  so¬ 
dium  vapor  based  on  resonance  saturation.  ,_5  Our  design 
differs  from  that  of  Boyd  etal.2  in  two  ways:  We  use  a  rectan¬ 


gular  frame,  which  permits  us  to  view  the  laser  produced 
plasma  channels  through  optically  flat  windows  and  we  have 
introduced  a  circular  thermal  expansion  buckle  to  relieve  the 
radial  stress  in  each  of  the  heated  plates.  A  cross-sectional 
view  of  one  half  of  this  heat  sandwich  oven  is  presented  as 
Fig.  1.  The  axis  of  symmetry  corresponds  to  the  extreme 
right  of  Fig.  1.  A  photograph  of  the  heat-sandwich  oven 
illuminated  by  the  ionizing  laser  pulse  is  shown  as  Fig.  2.  The 
narrow  sodium  plasma  channel  created  by  this  laser  is  clear¬ 
ly  visible  through  the  front  window. 

The  basic  frame  of  the  system  is  constructed  of  stainless 
steel  (type  304)  and  is  held  close  to  room  temperature  by 
circulating  water  through  the  cooling  coils  indicated  in  Fig. 
1 .  The  upper  and  lower  wicks  (required  for  heat  pipe  action) 
are  made  of  three  layers  of  stainless  steel  type  60  mesh 
(0. 1 52-mm  wires)  and  are  secured  to  the  heater  plates  (high- 
grade  stainless  steel-type  316)  by  a  large  number  of  small 
screws  and  to  the  cold  lip  by  a  ring.  The  rectangular  view¬ 
ports  were  each  20  cm  X  3  cm  in  size  and  the  heat  sandwich 
walls  had  O-ring  grooves  cut  into  them  in  order  to  provide  a 
seal  for  the  windows. 

Operation  is  commenced  by  pumping  down  the  argon 
pressure  until  it  is  just  below  the  sodium  partial  pressure 
required  for  a  given  experiment.  The  sodium  is  then  heated 
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Fit,  2  Photograph  of  sodium  plasma  created  by  laser  resonance  saturation 
within  the  heat  sandwich  oven 


until  it  melts,  whereupon  its  vapor  pressure  increases  and  the 
resulting  radial  flow  of  the  sodium  atoms  effectively  pumps 
the  argon  into  the  cool  outer  region  of  the  oven.  An  essential 
feature  of  this  heat  pipe  action  is  the  radially  inward  flow 
(from  cold  to  hot)  of  the  liquid  sodium  due  to  the  capillary 
forces  in  the  wicks.  It  should  be  noted  that  an  overpressure  of 
the  buffer  gas  leads  to  an  observed  drop  in  the  vapor  pressure 
of  the  sodium  for  under  these  conditions  the  argon  is  able  to 
permeate  the  hot  core  of  the  oven  and  thereby  transfer  ener¬ 
gy  from  the  hot  vapor  to  the  cool  boundary  region. 

Although  thermocouple  temperature  sensors  were  em¬ 
ployed,  it  was  felt  that  they  could  at  best  only  provide  a 
rough  estimate  of  the  sodium  atom  density  using  the  empiri¬ 
cal  fit  of  Nesmeyanov,h 


log,,,  X  =  29.84904  -  5619.4106/T 


-2.04111  iog,()r  +  3.45x10 -6r.  (1) 

As  one  can  see  from  Fig.  1,  the  thermocouple  wires  were 
placed  just  beneath  the  heaters.  Temperature  readings  from 
the  thermocouples  may  not  necessarily  reflect  that  of  the 
liquid  sodium  within  the  heat  sandwich  wick,  which  actually 
determines  the  vapor  pressure.  In  order  to  determine  the 
sodium  atom  distribution  across  the  oven  we  decided  to  em¬ 
ploy  the  absorption  measurement  technique  to  be  described 
below. 


II.  SODIUM  ATOM  DENSITY  MEASUREMENT 
TECHNIQUE 


The  spectral  irradiance  /(v,  l  |ofa  wcll-collimated  beam 
of  radiation  of  frequency  v  propagating  in  the  z-direction 
through  a  vapor  of  length  /  can  be  related  to  the  incident 
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spectral  irradiance  /  (v,  0)  by  the  Beer-Lambert  law . 

/(»-.  / 1  =  /  (v,  0)exp^  —  J~  k  |v.  zk/zj  ,  |2| 


where  k{v,z\  represents  the  volume  absorption  coefficient. 
In  the  event  that  v  is  close  to  the  resonance  doublet 
(3 :P  —  3 "S )  of  sodium  we  can  w  rite 


#r|v,  z)  --  — — .V  iJ”,  zi  Y  A  (J  ,  J  j  y  U  .  J ' 

Sir-  ■*— . 


V,  Zl  . 


I3l 


where  A  1/ ',  z)  represents  the  population  density  lat  a  depth  z 
into  the  vapor)  of  the  lower  level  having  a  total  angular  mo¬ 
mentum  quantum  number./ '  {J  '  —  1/2  for  sodium).  A  [J, .  J  j 
represents  the  Einstein  transition  probability  for  the  J,  to  J  ’ 
transition  of  the  doublet  (where  J-  =  3/2  and  J,  =  1/2  for 
the  resonance  doublet  of  sodiuml.  J  \ J, .  J  ,v,  z)  is  the  corre¬ 
sponding  line  profile  function,  and  /.  represents  the  mean 
wavelength  of  the  doublet. 

If  we  treat  the  absorption  as  arising  from  a  single  line 
with  an  "effective''  line  profile  function,  J  iv,  zi  we  can 
write 


8r 


V|/\  z)[A  (/,./  j  +  A  (/,./')]  J  iv. . 


where 

j  iv.  z)  =  :u,.j' \j  u,.J  '.v.  z) 

+  £(/;,  J'\J  (/;,  /',»',  Z)  ■  (5) 

In  this  instance  c(/,,/')and  ./(/,,/ ’,v.  z )  represent  the 
"fractional  line  strength"  and  “line  profile  function"  of  the 
J,  to/'  transition  of  the  doublet,  respectively. 

For  the  sake  of  brevity  we  introduce  the  simpler  nota¬ 
tion 

J  ,|v)  =  ._/  )/,,  J'.v.z) 

and 

A,  =  A  I /,,/')  for  /  =  1  or  2  . 

Also 

.V(z)  =  X(J  \  z\ . 

It  should  be  recognized  that  the  use  of  Eq.  (2)  implies  that  the 
radiation  is  well  below  the  saturated  value7  and  that  scatter¬ 
ing  is  negligible. 

At  the  600-800  K  temperature  range  corresponding  to 
the  1015  to  1017  cm  '  sodium  atom  density  within  the  heat 
sandwich  oven  both  Doppler  and  resonance  broadening 
have  to  be  considered.  Under  these  circumstances  the  resul¬ 
tant  resonance  D-line  absorption  profile  function  is  a  convo¬ 
lution  of  the  corresponding  Gaussian  |Doppler|  and  reso¬ 
nance  broadened  profiles,  that  is  to  say 


J 


•"-sM' 


dv *  e 


'  d  >  -  v*). 


(61 


where  v,  represents  the  ith  line  center  frequency,  J  fill  the 


resonance  broadened  atomic  line  profile,  and  fi,  the  Gaus¬ 
sian  width,  given  by  the  relation 

P,  =\2kT\-/mr\'  '- ■  (7! 

In  this  equation  k  and  c  have  their  usual  meaning.  T  is  the 
translational  temperature,  and  m  is  the  mass  of  the  sodium 
atom. 
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A.  Resonance  broadened  line  shape 


and 


The  core  of  the  resonance  broadened  atomic  line  profile, 
when  calculated  using  the  impact  approximation  (Appen¬ 
dix),  can  be  described  in  the  case  of  the  sodium  atom  by  a 
Lorentzian  with  half-width  at  half-maximum  (HWHMf 

T2  =  (1.47eVj/4jrv2m,)(Hz)  for(32S,,,  —  32/>,  ;) 
and  (8) 

T,  =  (1.805e2/,/4jrv,m,)(Hz)  for  (32S,/2  —  3 2P{  2I , 

where/  and  v,  are  the  oscillator  strength  and  line  center 
frequency  for  the  ;th  transition  and  me  and  e  being  the  mass 
and  charge  of  the  electron,  respectively. 

Srivastava  and  Zaidi0  have  shown  that  the  line  wings  of 
the  resonance  broadened  profile  can  also  be  described  by  a 
Lorentzian  with  HWHM  (when  calculated  using  the  quasi¬ 
static  approximation) 

T,  =  jT,  ,  (9) 

where  the  superscripts  w  and  c  denote  the  line  “wing”  and 
line  “core,"  respectively.  The  above  two  references  assume 
that  the  broadening  of  the  3 2P  state  is  dominated  by  the 
resonant  dipole-dipole  interaction  with  ground  state  sodium 
atoms. 

We  have  chosen  to  use  an  empirical  fit  to  the  line  profile 
of  the  form10 

.z»^c,f-i-ex  Pi-Iiv  — v^r.n 
Vr,  l  +  [(v-v,)/T,]2 


1  1  -exp]  -  [(v-v,.)/T,]2|  \ 

Vr,  1  +  j[(v  —  v,)/T,  ]2  )' 


which  approaches  the  exact  behavior  in  the  core  and  wings 
[v  =  i',  and  fv  —  v,  )/c  T,  >  1,  respectively]  and  resembles  the 
numerical  solution  in  the  intermediate  region.0  C,  is  the  nor¬ 
malization  constant 

C,  =  1.126/7,  (11) 

which  ensures 


Upon  substitution  of  Eq.  1 10)  into  Eq.  |6),  one  obtains  three 
integrals  which  all  involve  some  form  of  the  Voigt  integral 


W[  u,t]\  =  *  f  — - - -dy.  (13) 

tt-'tjJ  ,\^(u-yy 

We  have  chosen  to  use  the  approximation  of  Whiting"  for 


W(  U,  T) )  = 


.(■-irH 


4u:/w*  ] 


^11  ' 


J  W  ( u ,  y)du  =  1 .  (16) 

B.  Correction  for  finite  instrument  resolution 

In  our  experiments  a  Heath  scanning  monochromator 
was  employed  to  determine  the  spectral  characteristics  of 
the  transmitted  radiation.  For  this  kind  of  monochromator, 
where  the  entrance  and  exit  slits  are  comparable,  the  spectral 
transmission  function  T [vT  —  v)  is  well  approximated  by  a 
Guassian  profile  and  we  can  write 

7>r  -  V)  =  (1  /rz-T/'V'  -’-'V  (17) 

where  vT  is  the  frequency  of  peak  transmission  and 

yT  =5r/(ln2),/2.  (18) 

Here<5r  is  the  instrument’s  HWHM  bandwidth  correspond¬ 
ing  to  a  given  slit  setting  and  can  be  expressed  in  terms  of  the 
slope  bT  and  intercept  aT  of  the  linear  relation  between  the 
observed  HWHM  bandwidth  of  an  extremely  narrow  spec¬ 
tral  line  and  the  monochromator  slit  width  Sr, 

ST  =  aT  +  bTST  .  (19) 

The  observed  spectral  irradiance  at  frequency  vcan  thus 
be  written  in  the  form 


I(v*,l)T(v  —  v*)dv*  , 


or  using  Eqs.  (2)  and  (17)  and  assuming  broadband  incident 
radiation  so  that  7(v,0)=;/(0),  then  the  spectral  profile 

-  (v~  (2D 

n  ) 

can  be  solved  using  Eqs.  (4)— (6)  and  (10)  with  a  quadrature 
routine.  For  our  work,  the  slit  width  ST  was  150^m  and  the 
experimentally  determined  values  for  ar  and  bT  are  pro¬ 
vided  with  other  parameters  of  relevance  in  Table  I. 

To  evaluate  the  atom  density  distribution  it  is  necessary 
to  undertake  a  series  of  measurements.  The  sodium  vapor 
disk  is  assumed  to  be  divided  into  a  number  of  equally  spaced 
rings  of  constant  density  (see  Fig.  3)  and  the  absorption  is 
evaluated  for  a  series  of  chords  that  intercept  an  increasing 
number  of  these  rings.  Reference  to  Fig.  3  reveals  that  the 
element  of  path  length 

Tabi  r  I.  Parameters  of  relevance  in  sodium  atom  density  measurement 


4r/~  In  2l 


Vs 

=  5.09.1 

■:  10IJ  s  A 
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V, 

=  5.088 
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Fig.  3.  Schematic  diagram  illustrating  the  rings  of  constant  density  into 
which  the  sodium  vapor  disk  is  divided  for  the  purpose  of  determining  the 
radial  density  distribution. 


Ln  =  (r!+l  (ri  -X2m)l/2,  (22) 

where  xm  is  the  x  displacement  of  the  rath-measurement 
chord  and  r„  the  radius  of  the  nth-density  ring.  The  total 
optical  depth  for  the  rath-chord 

n  —  ^m«* 

r(v*,xj  =  2  2,  '  !23) 

n  =  1 

where 

*>*)  =  (A  VSrr )N{rn){A,  +  A2) (24) 

represents  the  absorption  coefficient  for  the  nth  ring  of  atom 
density  N \r„ )  and  (v*)  the  relevant  total  absorption  pro¬ 
file  for  the  sodium  D  doublet,  as  given  by  Eq.  (5).  In  order  to 
evaluate  the  Doppler  contribution  to  J?"  (v*)  for  each  ring 
the  temperature  is  assumed  to  be  related  to  the  atom  density 
N  ( rn )  of  that  ring  through  the  equilibrium  relation  given  ear¬ 
lier  as  Eq.  (1). 

The  density  for  the  nth  ring  is  then  guessed  and  the 
appropriate  absorption  spectrum  calculated;  this  is  then 
compared  with  the  corresponding  experimentally  deter¬ 
mined  spectrum.  The  assumed  density  is  then  adjusted  until 
reasonable  agreement  is  obtained.  The  average  density  of  the 
nth  ring  is  then  approximately  known  and  is  used  in  conjunc¬ 
tion  with  a  slightly  higher  density  for  the  two  segments  of  the 
next  ring  required  to  compute  the  absorption  spectrum  ex¬ 
pected  for  the  second  chord.  The  density  of  this  second  ring 
is  then  adjusted  until  reasonable  agreement  is  obtained 
between  the  calculated  spectrum  and  that  measured.  This 
series  of  calculations  and  comparisons  with  the  measure¬ 
ments  is  continued  towards  the  center  of  the  heat  sandwich 
oven  leading  to  an  effective  radial  density  profile  for  the  sodi¬ 
um  vapor. 

Reasonable  agreement  is  obtained  when  the  full-width 


f 

I 

I 

at  half-maximum  (FWHM)  of  the  predicted  absorption  spec¬ 
trum  closely  matches  that  measured  well  within  the  limit  of 
instrument  resolution  (see  top  inset  of  Fig.  5).  When  the  sodi¬ 
um  doublet  is  clearly  resolvable  (usually  for  x  £  5  cm),  we 
used  the  absorption  at  line  center  of  the  32P3/2  —  32S1/2  tran¬ 
sition  as  our  matching  condition  (see  lower  inset  of  Fig.  5).  It  j 
should  be  noted  that  best  results  were  obtained  when  pre¬ 
liminary  smoothing  of  the  data  was  performed.  The  FWHM 
of  the  absorption  spectrum  was  fit  to  a  polynomial  in  posi¬ 
tion  x  of  the  form  , 

P(x)  =  (a,  +  a^x2  +  a3jc4)  exp  (  -  a4x 2) ,  (25) 

where  nr,  (/'  =  1  -4)  represents  the  various  fitting  parameters.  ! 
The  variation  of  the  transmission  at  line  center  of  the  D2  line  ' 

with  position  x  was  approximated  by  a  straight  line.  Data 
from  these  two  curves  were  then  directly  used  in  the  density 

reconstruction  routine.  1 

i 

III.  EXPERIMENTAL  ARRANGEMENT  AND  RESULTS 

The  facility  for  measuring  the  sodium  atom  distribution 
comprised  a  broadband  xenon  light  source,  an  aperture  and 
collimating  optics,  the  heat  sandwich  oven,  a  Heath  model 
EU-700  scanning  grating  (1180  lines/mm,  / /6.8,  35-cm  fo¬ 
cal  length)  monochromator,  and  an  RCA  type  1P28  photo¬ 
multiplier.  A  schematic  of  the  optical  arrangement  is  pre¬ 
sented  as  Fig.  4.  The  diameter  of  the  probe  beam  of 
continuum  radiation  was  set  to  about  0.3  cm  by  the  stop  S. 

An  example  of  the  variation  of  the  FWHM  and  D2  line 
center  transmission  with  position  x  is  presented  as  Fig.  5. 

The  two  insets  display  two  representative  absorption  pro- 


Fig  4  Facility  for  determining  the  sodium  vapor  density  by  means  of  mea¬ 
suring  the  absorption  spectrum  The  lenses  arc  designated  L.  aperture  S, 
pinhole  PH.  and  mirrors  M 
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Pressure  ;  36  mm 


40  6C 


Fig.  5  Vanation  of  the  full-width  half-maximum  Isolid  linel  and  D.  line 
center  transmission  (dashed  hnel  with  position  X.  In  this  run  the  argon  pres¬ 
sure  was  about  3.6  mm  of  mercury  and  the  maximum  temperature  was  938 
K.  The  two  insets  display  examples  of  the  absorption  spectrum  measured 
for  a  chord  passing  through  the  core  of  the  vapor  disk  and  one  passing 
through  the  low -density  rim. 


files,  one  corresponding  to  a  path  close  to  the  center  of  the 
oven  the  other  in  the  low-density  rim  of  the  vapor  disk.  In  the 
latter  the  absorption  arising  from  the  two  sodium  D  lines  is 
evident,  while  in  the  former  the  optical  depth  is  so  great  that 
only  one  wide  spectral  hole  is  apparent.  The  corresponding 
radial  density  distribution  of  sodium  atoms  within  the  heat 
sandwich  oven  is  presented  as  Fig.  6.  This  distribution  was 
determined  from  a  series  of  absorption  measurements  in  ac¬ 
cordance  with  the  technique  described  earlier.  The  maxi¬ 
mum  temperature  registered  by  the  thermocouple  near  the 
oven  center  was  938  K,  which  would  suggest  a  peak  density 
considerably  larger  than  that  measured.  Mechanical  equilib¬ 
rium  between  the  argon  buffer  gas  and  sodium  vapor  was 
ensured  by  continuously  monitoring  and  adjusting  the  buff- 


Tmax  -  93F  K 
Pressure  -  36  mm 


Fig.  6  Sodium  atom  density  distribution  within  the  heat  sandwich  oven 
(corresponding  to  the  data  of  Fig.  5|  determined  by  absorption  and  the  dc 
convolution  technique  described  in  the  text 


r  ;  cm  ) 


Fig.  7.  The  sodium  atom  distribution  measured  within  the  heat  sandwich 
oven  for  an  argon  pressure  of  3.05  mm  of  mercury  and  a  maximum  tempera¬ 
ture  of  about  903  K 

er  gas  pressure  to  match  the  sodium  vapor  pressure  at  the 
center  of  the  heat  sandwich  oven. 

An  important  observation  to  be  drawn  from  this  density 
profile  is  that  the  radial  atom  distribution  proceeds  to  fall  by 
an  order  of  magnitude  over  the  first  4  cm  of  radius.  A  similar 
result  was  obtained  when  the  oven  was  operated  at  a  lower 
temperature  such  that  the  peak  sodium  density  was  about  a 
factor  of  3  less  than  indicated  for  the  case  shown  in  Fig.  6. 
The  atom  density  distribution  in  this  instance  is  presented  in 
Fig.  7. 


0  2  9  6  8 


r  cm 

Fig.  8  The  sodium  atom  distribution  measured  within  the  heat  sandwich 
oven  for  an  argon  pressure  of  4  0  mm  of  mercury  and  a  maximum  tempera¬ 
ture  of  about  953  K. 
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These  results  are  to  some  extent  in  keeping  with  the 
analysis  of  Vidal12  and  Boyd  and  Dodd11  who  suggests  that 
substantial  variations  of  density  are  possible  in  heat-pipe  sys¬ 
tems  as  a  result  of  the  large  flow  velocities  that  can  be  estab¬ 
lished. 

We  have  found  that  operation  of  the  heat  sandwich  oven 
for  an  extensive  period  of  time  |  £  10  h)  at  relatively  high 
temperatures,  with  an  underpressure  of  argon  gas,  leads  to  a 
“donut”-shaped  distribution  of  sodium  vapor.  Figure  8  dis¬ 
plays  such  a  distribution.  We  believe  that  this  arises  from  too 
great  a  heater  power  combined  with  too  low  a  buffer  gas 
pressure.  Under  these  conditions,  large  radial  flow  velocities 
were  established  and  the  sodium  was  found  to  condense  and 
solidify  beyond  the  region  of  the  wick,  resulting  in  a  loss  of 
sodium  and  a  gradual  reduction  in  the  amount  available  to 
fill  the  heated  volume.  In  time,  the  center  of  the  wick  dries 
out  and  any  liquid  sodium  which  migrates  inwards  from  the 
condensation  zone  as  a  result  of  the  capillary  action  would 
evaporate  before  reaching  the  oven  center.  We  assume  that 
under  these  circumstances  a  pressure  balance  is  maintained 
by  the  formation  of  an  argon  bubble  in  the  center  of  the  heat 
sandwich  oven. 

These  surprising  results  suggest  that  operation  of  the 
oven  in  a  heat  pipe  mode  requires  a  continuous  monitoring 
of  both  the  vapor  density  distribution  and  buffer  gas  pres¬ 
sure.  Periodic  adjustment  of  the  buffer  gas  pressure  may  be 
necessary  to  avoid  the  creation  of  this  donut-like  distribu¬ 
tion. 

We  estimate  the  accuracy  of  the  atom  density  measure¬ 
ments  presented  here  to  be  within  ±  20%.  This  is  primarily 
determined  by  the  uncertainty  in  the  half-width  of  the  reso¬ 
nance  broadened  atomic  absorption  profile.  The  theoretical 
values  by  Carrington  et  al.H  used  here  agree  fairly  well  with 
the  experimental  measurements  of  Huennekens  and  Gal¬ 
lagher14  which  were  estimated  to  have  an  error  of  no  more 
than  ±  15%.  We  have  found  that  in  the  density  regions  of 
interest,  there  is  at  most  a  direct  relationship  between  the 
percentage  change  in  the  half-width  of  the  atomic  absorption 
profile  and  the  percentage  change  in  the  peak  atom  density 
evaluated. 

Secondary  sources  of  error  are  contributions  to  the 
atomic  line  profile  due  to  van  der  Waals  interactions  (Ap¬ 
pendix),  as  well  as  fluctuations  in  the  oven  temperature  (  ±  5 
K  leading  to  ;=  10%  density  variation)  during  the  series  of 
measurements.  There  was  no  noticeable  change  in  the  out¬ 
put  of  the  xenon  lamp  between  each  scan.  Currently  we  are 
attempting  to  reduce  the  uncertainty  as  a  result  of  tempera¬ 
ture  fluctuations  by  making  these  density  measurements  in  a 
much  shorter  time. 
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APPENDIX 

If  absorption  of  a  photon  by  an  atom  takes  place  during 
a  collision  of  that  atom  with  another  atom  Iperturber),  then 


the  transition  can  be  thought  of  as  occuring  between  two 
quasimolecular  states  formed  as  a  result  of  the  interaction 
between  the  wavefunctions  of  the  colliding  partners.  The  re¬ 
sulting  energy  shifts  of  the  excited  and  ground  states  are  a 
function  of  the  interatomic  separation  and  depend  on  the 
type  of  interaction.  The  curves  depicting  the  potential  ener¬ 
gy  of  these  quasimolecular  states  as  a  function  of  intemu- 
clear  separation  r  are  referred  to  as  the  ground  or  excited 
state  potentials. 

For  the  case  of  sodium  perturbed  by  another  sodium 
atom,  the  3  2P  states  of  the  absorbing  atom  split  up  into  four 
potential  curves1'  (neglecting  spin-orbit  and  spin-spin  inter¬ 
actions)  situated  symmetrically  about  the  atomic  3 2P  state, 
at  large  intemuclear  separations  (long  range).  At  long  range, 
the  splitting  of  the  3 2P  state  primarily  determines  the  shape 
of  the  atomic  absorption  profile  and  is  dominated  by  the 
long-range  resonant  dipole-dipole  interaction. 16  The  magni¬ 
tude  of  the  energy  shift  has  the  form 

JtAco  =  C,/r' ,  (Al) 

where  C}/h  is  the  resonant  dipole-dipole  interaction  con¬ 
stant,  A co  =  (o  —  (o0,  co  the  frequency  of  the  absorbed  pho¬ 
ton,  co0  that  in  the  absence  of  a  perturber,  and  r  is  the  intemu¬ 
clear  separation.  At  very  large  intemuclear  separations, 
where  the  duration  of  the  collision 


Td  =Rc/v  (A2) 

(here  Rc  is  the  impact  parameter  and  u  the  relative  velocity) 
is  much  less  than  the  correlation  time  describing  the  phase 
relationship  between  two  points  in  time  of  the  amplitude  of 
the  electron  treated  as  a  classical  oscillator,17  i.e., 

rrf  <l/(<y  —  &>,,),  (A3) 

then  the  effect  of  the  collision  can  be  treated  as  a  sudden 
disruption  of  phase  in  the  amplitude  of  the  oscillator.  Treat¬ 
ed  under  this  “impact”  approximation,  the  absorption  pro¬ 
file  is  Lorentzian  in  shape 


s’  M  = 


2n  («  —  cog)2  +  y2/ 4 


with  full-width  half-maximum  (FWHM)  y  varying  linearly 
with  perturber  density  N.  The  value  y/N  is  often  referred  to 
as  the  resonance  broadening  parameter.1' 

In  the  other  extreme,  where 


Td>\/(co  —  <u„),  [A5| 

the  field  experienced  by  the  absorber  as  a  result  of  a  neigh¬ 
boring  perturber  a  distance  r,  away  is  assumed  constant  in 
time,  and  the  line  shape  is  proportional  to  the  probability  of 
finding  the  nearest  particle  a  distance  r  from  the  absorber. 
Assuming  that  the  nearest  perturber  is  much  closer  than  the 
average  separation  between  atoms  (nearest  neighbor  or  bina¬ 
ry  approximation),  yet  sufficiently  far  apart  such  that  the 
interaction  does  not  significantly  alter  the  distribution  of 
atoms  about  the  absorber,  then  we  can  describe  the  probabil¬ 
ity  P[r\dr  of  the  nearest  perturber  being  within  a  distance  r 
and  r  +  dr  from  the  absorbing  atom1'’ 

P(r]dr  =  4 vrS  exp[  -  (4n-/3)Ar']<fr  .  |A6| 

According  to  the  statistical  theory,  the  line  shape 

J  \it)\do>  ■’  P[r\dr  i  A7) 
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and  substituting  (Al)  into  (A6),  we  arrive  at 


$ 

£4 


jf'{a))dco  <x 


4  nNCy 


-  ■> « 
V> 


3Aor  rV  A  a>F  ' 

where  we  have  defined 

hAa0  =  Cy/r^  (A9) 

with  r0  representing  the  mean  intemuclear  separation,  i.e., 

r0  =  (3/4irJV )'■'•' .  (A  10) 

Applying  the  conditions  that 

j  J\(o\da>  =  (  (All) 

Jo 

[since  the  splitting  is  in  the  form  of  ±  C,/r*  (Ref.  15)  and 
that  we  are  only  considering  one  side  of  the  line  profile]  and 
assume  that  A<u0<A io  (or  r<r(1,  i.e.,  nearest-neighbor  ap¬ 
proximation),  then 

XM  =  |2jr/3|.VC,(Aft»r-’  (A12) 

and  the  wings  of  the  profile  are  also  Lorentzian  in  shape. 

Of  course,  this  “quasistatic"  approximation  is  inappli¬ 
cable  to  the  line  core  as  one  can  see  that  it  would  violate  the 
conditions  described  by  (A5).  In  the  line  core,  the  shape  is 
described  by  Eq.  (A4)  which  conversely  is  inapplicable  to  the 
line  wings. 

Detailed  calculations  of  the  resonance  broadened  line 
profile  describing  the  line  wings,  core,  and  the  intermediate 
regions  where  rd  ~  \/(oj  —  oQ) ,  have  been  performed  by  Sri- 
vastava  and  Zaidi.9  They  have  found  that  the  effective 
FWHM  of  the  Lorentzian  describing  the  quasistatic  wings  is 
two-thirds  that  of  the  profile  describing  the  line  core 

wr  =  Vr>  (Ai3) 

where  the  superscripts  c  and  w  refer  to  the  core  and  wings, 
respectively. 

Under  the  range  of  conditions  expected  for  our  experi¬ 
ments,  with  maximum  sodium  density  Amax  =  5  X 1016 
cm  " 3,  the  average  intemuclear  separation 

r0  =  (3/47rAmaJ,/2=;168  A.  (A14) 

The  binary  approximation  describing  the  distribution  of 
the  nearest  perturber  is  valid  for  r<r0.  At  an  interatomic 
separation  of  /■</ 10,  it  remains  to  check  the  frequency  range 
imposed  by  the  criteria  for  applicability  of  the  quasistatic 
approximation  [Eq.  (A5)] 

A<y>10i ;/r0.  (A  15) 

For  relative  velocities  of  approximately  105  cm/s,  we  arrive 
at 

A<u>7x  10"  rad/s  ,  {A  1 6) 

which,  for  the  sodium  resonance  lines,  corresponds  to 

AZ>1  A  .  (A17) 

Generally,  when  passing  the  broadband  radiation  through 
the  oven  center  spectral  holes  of  greater  than  25  A  were 
observed.  For  such  measurements,  the  sodium  density  is  es¬ 
timated  from  the  absorption  at  approximately  1 0  A  from  line 
center,  thereby  satisfying  the  criterion  imposed  by  (A17). 

Where  the  broadband  radiation  passes  through  regions 
of  low  density  and  short  optical  path  (i.e.,  an  outer  chord  of 
the  vapor  disk),  the  density  is  estimated  from  the  absorption 


at  line  center  of  the  32A,/2  —  3JP3/2  transition.  The  line  cen¬ 
ter  of  course  satisfies  the  criteria  of  (A3),  thereby  permitting 
the  use  of  the  impact  approximation  to  describe  the  shape  of 
the  line  core. 

At  shorter  interatomic  separations  ( ~r(/10)  one  should 
be  cautious  to  check  that  the  van  der  Waals  dipole-dipole 
interaction  plays  a  negligible  role  in  the  perturbation  of  the 
3  2P  state.  This  interaction  leads  to  an  energy  shift  of  the  state 

hkco  =  Cb/rb.  (A  18) 

We  have  estimated  the  C6  contribution  to  the  excited  stated 
'  2  *  van  der  W aals  potential  by  fitting  the  numerically  tabu¬ 
lated  potentials  by  Kaminsky18  to 

AA*=  (A19) 

r  r  r 

where  we  have  also  included  the  quadrupole-quadrupole 
CgZ-8  interaction  term.  Using  the  value  of  C,  experimentally 
determined  by  Niemax  and  Pichler  l5c 

Cyh  =  5  X  1CT 8  cm3  s  ~ 1 ,  (A20) 

we  arrive  at 

Cf/h  =  5x  10“ 29  cm"  s- ',  Cg/h  =  10-43  cm8  s_l . 

(A21) 

The  shift  of  energy  state  due  to  the  resonance  and  van 
der  Waals  interaction  becomes  comparable  at  interatomic 
separations 

r  =  (C6/C3),n  =  10  A  .  (A22) 

Using  Eqs.  (A19HA21)  with  r  =  10~7  cm  (10  A),  then 
(a?  —  co0)  =  1.1  X  1014  rad  s_1  (A23) 

or 

(A  -  A0)  =  202  A,  (A24) 

at  which  one  can  no  longer  exclude  the  van  der  Waals  inter¬ 
action  in  determining  the  shape  of  the  extreme  line  wings.  If 
we  limit  our  measurements  to  approximately  ±  20  A 
(q)—o) 0=1.1x10' 3  rads-1)  from  either  resonance  line 
center,  then  Eq.  (A19)  leads  to 

r~  18  A  (A25) 

and  at  such  intemuclear  separations,  the  contribution  to  the 
potential  due  to  the  van  der  Waals  interaction  is  about  15% 
that  due  to  the  resonant  dipole-dipole  interaction.  This 
would  translate  to  a  15%  change  in  the  absorption  coeffi¬ 
cient  which  we  have  found  to  alter  the  peak  sodium  density 
measured  by  approximately  the  same  amount.  Only  in  few 
cases  however  have  we  operated  the  oven  at  such  high  densi¬ 
ties  and  have  measured  such  large  spectral  holes. 

van  der  Waals  broadening  due  to  a  foreign  gas 

At  some  intemuclear  separation  r,  the  contribution  to 
the  potential  due  to  the  van  der  Waals  interaction  with  a 
foreign  gas  as  a  perturber  would  be 

AA«  =  C{/r6,  (A26) 

where  C{/h  is  the  van  der  Waals  interaction  constant  esti¬ 
mated  by  Keilkopf'9  for  argon,  to  be 

C{/h~5x  10-30cmAs_l ,  (A27) 
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which,  for  intemuclear  separations  of  approximately  1 8  A, 
leads  to 

[co  —  eo0)  =  1.5  X  10"  rad  s~‘ ,  (A28) 

which  is  orders  of  magnitude  less  than  that  resulting  from 
the  resonant  dipole-dipole  interaction.  At  argon  pressures 
comparable  to  the  sodium  vapor  pressure,  we  could  safely 
ignore  the  broadening  of  the  atomic  absorption  profile  due  to 
argon  present  as  a  buffer  gas. 
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Abstract 

The  first  free  electron  density  radial  profiles  of  a  sodium  plasma 
created  by  laser  resonance  saturation  are  reported.  The  measurements  were 
based  on  Stark  broadening  of  the  42D-32P  multiplet  and  reveal  that  strong 
absorption  of  the  laser  pulse  leads  to  formation  of  a  conically  shaped 
plasma  along  the  path  of  the  laser  pulse. 


INTRODUCTION 

Electron  density  measurements  from  the  Stark  broadening  of  spectral 

emission  lines  has  developed  into  a  useful  and  common  diagnostic  tool  since 

1-4 

the  formulation  of  reliable  Stark  broadening  theories.  The  most 
attractive  feature  of  this  technique  is  that  it  is  noninvasive  and  has  been 
widely  used  in  applications  where  plasma  ionization  probes  either  directly 
interfere  with  the  processes  studied,  or  are  easiy  contaminated  by  highly 
reactive  constituents  of  the  plasma.  The  viability  of  the  technique  was 

5 

clearly  demonstrated  by  Agnew  and  Reichelt  in  determining  the  electron 

density  in  a  cesium  plasma  diode,  and  afterwards  in  wall  stabilized  arcs6"10 

where  cylindrical  symmetry  and  plasma  stability  permitted  the  Abel  inversion 

of  lateral  intensity  measurements  to  obtain  the  local  emission  coefficient, 

and  hence  the  radial  variation  in  electron  density.  Concurrently,  Stark 

broadening  has  been  widely  investigated  and  employed  as  a  diagnostic  tool  in 

11-14 

transient  plasmas  produced  in  shock  tubes,  laser  ablation  of  solid 

targets,16,16  low  induction  vacuum  sparks1^’10  and  in  plasmas  produced  by 

19-21 

laser  resonance  saturation.  Most  of  these,  however,  were  spatially 

and/or  temporally  averaged,  depending  on  the  nature  or  reproducibility  of 

20  21 

the  source.  In  some  cases,  multi  shot  averaging  has  been  performed.  * 


beam.  The  maynif Ration  of  the  optics  was  1:1  with  the  SPEX  entrance  slit 


(Academic,  New  'Torn,  1962). 
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Abstract 

A  computational  study  of  the  effects  of:  inhomogeneities, 

self-absorption  and  finite  detection  bandwidth  on  the  electron  temperatures 
deduced  from  Boltzmann  plots  of  line  intensities  from  sodium  plasmas  in  LTE 
has  been  undertaken  and  the  results  compared  with  an  experiment, 

INTRODUCTION 

The  measurements  of  the  electron  temperatures  from  the  relative 

1-4 

intensities  of  spectral  lines  is  a  common  approach  in  plasmas  where 
either  an  absolute  measurement  of  line  intensity  is  not  possible,  or  the 
recombination  radiation  is  weak  and  where  limited  spectral  resolution 
prevents  the  measurement  of  the  shift  to  width  ratio  of  Stark  broadened 
emission  lines.  The  applicability  of  a  “Boltzmann  plot",  the  slope  of 
which  can  be  used  to  evaluate  the  free  electron  temperature,  T  ,  in  an  LTE 
plasma  and  the  simplicity  of  the  technique,  rests  on  the  assumptions  that 
the  plasma  is  homogeneous  and  that  the  spectral  lines  suffer  negligible 
sel  f-absorption. 

The  first  assumption  is  rarely  justified  and  even  in  the  case  of 
cylindrical  symmetry,  a  more  detailed  analysis  is  required,  involving  the 

g 

solution  of  the  Abel  transform  to  reconstruct  the  radial  volume  emission 

g 

coefficients  from  lateral  intensity  measurements.  The  latter  assumption  is 
generally  valild  for  lines  with  high  lying  upper  states  whose  terminating 
state  is  any  state  other  than  the  ground  state.  Such  lines  have  relatively 
weak  oscillator  strengths  and  are  usually  optically  thin.  These  same  lines. 


slopes  and  therefore  the  predicted  temperatures  differ  somewhat  between  the 
two  configurations  (by  rouyhly  10%).  Where  the  temperature  yradients  are 
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Computed  Boltzmann  Plots  for  Several  Transitions  in  a  Sodium 
Plasma  with  Flat  Radial  Distributions  of  Electron  Temperature 
Corresponding  to:  (D  lO.OOOK  and  (II)  5.000K  (Shown  in  the 
Small  Inset). 
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Radial  Distributions  of  40  Population  in  a  Sodium  Plasma 
Corresponding  to  the  Two  Radial  Electron  Temperature  Profiles 
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tzmann  Plots  for  Several  Transitions  Expected  from  Conf iguration 
or  a  Sodium  Plasma  Assumed  to  have  the  Two  Temperature 
tributions  Indicated  in  Figure  4  (a)  and  (b).  Both  Optically 
n  and  Full  Solution  Data  are  Presented.  Solutions  Corresponding 
Observation  with  a  3nm  Detection  Bandwidth  are  also  Presented. 
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Abstract 

The  suitability  of  assuming  a  Gaussian  radial  profile  for  the  incident 
laser  pulse  used  in  a  computational  model  of  a  highly  nonlinear  interaction 
(laser  resonance  ionization  of  sodium  vapor)  is  addressed  when  the 
predictions  of  this  model  are  to  be  compared  with  the  multishot  average 
results  of  an  experiment  using  radially  noisy  incident  laser  pulses. 

INTRODUCTION 

The  experimental  measurements  of  spatial  distributions  of  quantities 
such  as  electron  density  and  temperature  in  a  plasma  are  frequently  derived 
from  an  analysis  of  the  spatially  integrated  spectral  radiances  arising  from 
high  lying  atomic  or  ionic  transitions.  In  many  pulsed  experiments  the 
spectral  distribution  of  the  emission  is  determined  from  a  series  of 
measurements  at  different  wavelengths.  If  the  shot-to-shot  reproducibility 
is  poor,  or  the  signal  is  weak,  the  emission  recorded  in  each  spectral 
interval  is  an  average  of  several  shots.  When  the  results  of  such 
experiments  are  to  be  compared  with  those  predicted  by  a  theoretical  model, 
difficulties  can  arise  if  the  shot-to-shot  variation  of  the  interaction 
responsible  for  the  observed  emission  is  appreciable  and  the  dependence  is 
highly  nonlinear. 

This  situation  is  well  illustrated  in  the  case  of  plasmas  created  by 

laser  resonance  saturation.  The  ionization  capability  of  resonantly  tuned 

(1  2) 

laser  radiation  was  experimentally  discovered  by  Lucatorto  and  Mcllrath  ’ 

(3  5) 

and  theoretically  accounted  for  by  Measures  and  his  coworkers.'  ’ 
However,  the  strong  attenuation  of  the  laser  pulse  and  the  highly  nonlinear 
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Figure  4.  Computer  Predicted  Redial  Variation  of  (a)  Free  electron  Density 
and  Temperature,  (b)  4D  Population  Density  at  t=58ns  and  Z=0. 
Also  shown  is  the  Noisy  Radial  Profile  of  the  Incident  Laser 
Pulse  Used  in  the  Computation. 
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ABSTRACT 


We  have  developed  a  computer  cuor  that  models  the  three  dimensional 
nature  of  a  sodium  plasma  created  try  laser  resonance  saturation.  This  was 
accomplished  by  taking  account  of  the  temporal  and  spatial  distortion 
suffered  by  the  laser  pulse  as  it  propagates  through  a  nonuniform  sodium 
atom  distribution.  We  have  also  undertaken  the  first  measurements  of  the 
radial  and  axial  electron  density  and  temperature  profiles  arising  in  such 
a  plasma  and  have  been  able  to  demonstrate  that  our  computer  code 
predictions  are  in  reasonable  agreement  with  these  results.  Absorption  of 
the  laser  ulse  is  found  to  be  quite  significant  for  even  a  few  centimeters 
of  1016  cm-3  sodium  atom  density  and  that  as  a  result  quite  low  temperature 
plasmas  can  be  produced  for  modest  energy  laser  pulses. 


INTRODUCTION 

Last"  saturation  of  an  atomic  resonance  transition  clearly  involves 

radiative!,)  locking  the  resonance  and  ground  state  populations  in  the  ratio 

(11  (2  31 

of  their  degeneracies.  '  Lucatorto  and  Mcllrath'  *  '  were  the  first  to 

show  experimental ly  that  this  interaction  represents  an  effective  method  of 

producing  near  total  ionization  of  alkali  atoms.  Subsequent  experiments 

have  not  only  verified  this  observation  but  have  shown  that  this  form  of 

(4-71 

laser  ionization  can  be  extended  to  other  elements/  '  There  is  growing 

evidence  that  absorption  and  superelastic  electron  quenching  collisions 

together  represent  the  primary  mechanism  for  rapidly  extracting  from  the 

laser  pulse  the  large  amount  of  energy  needed  to  achieve  this  hiyh  degree 

of  ionization.  This  mechanism  was  first  suggested  by  Measures, 

and  is  schematically  illustrated  in  figure  1.  A  useful  overview  of  the 

(51 

subject  was  provided  by  Jahreiss  and  Huber/  ' 

( 19-21 1 

Measures  and  coworkers/  have  proposed  that  laser  ionization 

based  on  resonance  saturation  (LIBORS)  proceeds  in  four  stages:  (i)  the 
laser  rapidly  "locks"  the  ground  and  resonance  level  populations  roughly  in 
the  ratio  of  their  respective  degeneracies,  (ii)  there  is  a  rapid  yrowth  of 
free  electrons  due  to  one  or  more  "seed  ionization"  processes  and  these 
electrons  very  quickly  acquire  energy  through  "superelastic  collision 
quenching"  of  the  laser  maintained  resonance  state  population,  (iii)  these 
"hot"  electrons  proceed  to  col  1 isional ly  excite  and  ionize  the  resonance 
state  atoms  and  in  addition  the  col  1 i sional ly  excited  higher  lying  levels 
are  "photoi oni zed"  by  the  laser  field,  finally  (iv)  "runaway"  collisional 
ionization  of  the  upper  levels  occurs  once  some  critical  free  electron 
density  is  achieved  and  this  leads  to  near  complete  ionization  burnout  of 
the  species. 

■  However,  it  should  be  noted  that  if  a  very  high  rate  of  seed 

ionization  exists  the  4  states  telescope  to  3  and  the  burnout  stage 

(191 

precedes  almost  directly  from  seed  ionization/  ' 

During  the  ionization  burnout  stage  superelastic  heating  can  no  longer 

balance  collisional  cooling  due  to  depletion  of  the  resonance  state 

( 191 

population  and  consequently  the  electron  temperature  is  predicted'  '  to 
fall  significantly.  This  effect  has  also  been  predicted  by  Morgan^22) 
who  has  modelled  the  ionization  using  the  Boltzmann  equation  rather  than 
assuming  that  the  free  electrons  have  a  Maxwellian  velocity  distribution. 


In  the  original  formulation  of  our  LIBORS  collisional -radiative 
computer  code  the  sodium  atom  was  represented  k  c  2u-atomic  level  model 

( 1  q  _  p  i  ;>3  i 

and  we  assumed  a  step-like  laser  pulse.  -<-»-<  An  example  of  the 

temporal  behaviour  of  the  free  electron  density,  free  electron  temperature, 

resonance  state  density,  and  ground  state  density  predicted  by  this  code 

modified  to  include  a  realistic  flashlamp  pumped  dye  laser  pulse  (such  as 

used  in  the  early  experiments  ^ )is  presented  in  figure  2.  Also  shown  on 

this  figure  is  the  assumed  temporal  growth  of  the  laser  irradiance.  The 

influence  of  molecular  nitrogen  on  this  ionization  process  has  also  been 

(241  (251 

considered  by  Measures  et  al '  and  Wong.'  '  In  an  actual  experiment  the 
laser  pulse  propagates  through  a  non-uniform  atom  density  distribution  and 
suffers  appreciable  absorption  due  to  the  strong  nature  of  the  interaction. 
This  can  lead  to  a  considerable  temporal  and  spatial  distortion  of  the 
laser  pulse.  Consequently,  the  atoms  at  different  locations  along  the  path 
of  the  laser  pulse  are  subject  to  progressively  less  energy  and  to  a 

shorter  period  of  excitation. 

In  order  to  better  model  the  real  effects  associated  with  creating  a 
plasma  channel  by  laser  resonance  saturation  we  have  undertaken  the  first 
analysis  that  takes  account  of  both  laser  absorption  and  the  non-uniform 
atom  density  distribution  along  the  path  of  the  laser  pulse.  We  have  also 
taken  into  account  the  spatial  and  temporal  profile  of  the  laser  pulse. 

In  order  to  test  this  theoretical  work  we  have  attempted  to 
computationally  simulate  as  closely  as  possible  experiments  that  we  have 


undertaken  using  a  specially  designed  sodium  heat  sandwich  oven 
Nd-YAG  laser  pumped  dye  laser.  x 


(26) 


and  a 


5-LEVEL  ATOM  MODEL 


In  an  attempt  to  keep  the  analysis  economically  viable  we  replaced  our 

(19-21  23) 

20-level  model  of  the  sodium  atom'  *  '  with  the  simplified  5-level 

( 25 ) 

model'  '  shown  in  figure  3.  Although  5  levels  are  indicated  in  figure  3, 
we  assume  that  the  collisional  coupling  between  levels  3  (42S),  4  (32D)  and 
5  (42P)  is  sufficiently  rapid  that  their  populations  are  maintained  in  a 
Boltzmann  equilibrium.  This  assumption  permits  us  to  treat  this  trio  of 
levels  as  a  manifold  and  thereby  reduce  the  number  of  atomic  level  rate 


E.4 


equations  to  three.  In  effect  we  consider  only  the  collisional-radiative 
processes  "into"  and  "out  of"  tnis  (3-4-5)  level  manifold. 

The  remaining  15  high  lying  energy  levels  used  in  tne  ?j- level  model 
(not  shown  in  figure  3)  are  partially  taken  into  account  by  assuming  that 
electron  excitation  to  them  (from  any  of  the  five  lower  levels)  leads,  in 
effect,  to  immediate  ionization  of  the  atom.  We  therefore  introduce  the 
5-level  "effective"  electron  collisional  ionization  rate  coefficient  for 
each  of  the  low  lying  levels  (m  =  1  to  5)  defined  by  the  relation 

■Sc  ■  X  Km  (D 

n=o 

where  K  represents  the  electron  collisional  rate  coefficient  between 
mn  K 

levels  m  and  n.  Note  that  K^jis  the  m_1evel  electron  collisional 
ionization  rate  coefficient  according  to  the  20-level  model  of  the  sodium 
atom.  The  corresponding  "effective"  three  body  recombination  rate 
coefficient  to  each  of  the  low  lying  levels  (m  =  1  to  5)  is  similarly 
constituted  to  take  account  of  recombination  to  the  higher  levels,  followed 
by  a  rapid  cascade  down  to  the  lower  levels.  We  thus  introduce 

20 

Kc*m  *  \  KcnP(n-ra)  (Z> 

n=o 

where  Kcn  represents  the  three-body  recombination  Tate  coefficient  to  level 

n  used  in  our  20-level  LIBORS  code  and  p(n,m)  is  introduced  as  the 

* 

effective  probability  that  a  captured  free  electron  arriving  in  level  n 

subsequently  decays  to  level  m.  In  principle  radiative  recombination  could 

be  included  into  the  effective  recombination  rates  but  for  our  conditions 

(25) 

its  omission  was  found  to  lead  to  negligible  error. 

In  calculating  the  electronic  excitation  rate  coefficients  we  have 

(271 

averaged  Seaton  cross-sections'  for  the  optically  allowed  transitions, 

and  Gryzinski  cross-sections^*^  for  all  other  transitions,  over  a 

Maxwellian  velocity  distribution.  In  the  case  of  the  resonance  transition) 

( 29) 

the  experimental  value  of  Crandall  et  al ,  was  employed.  Although  both 
Seaton  and  Gryzinski  collision  rates  were  employed  for  transitions  between 


the  5  lower  levels  only  Seaton  collision  rates  were  used  for  transitions 
that  terminated  above  level  5. 


We  have  found  that  close  agreement  can  be  obtained  between  the  5  and 
20-level  LIBORS  codes  if  these  probability  factors  are  generallyt  given  by 
the  expression 


P(n,m) 
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where  N  is  the  free  electron  density,  A  the  Einstein  spontaneous 
e  J  nm  r 

transition  probability  between  levels  n  and  m,  with  q  being  the  levels 
within  the  manifold  that  are  optically  connected  to  level  n. 

The  5-level  model  can  thus  be  described  by  the  following  set  of 
equations:  first  we  have  the  three  atomic  level  population  rate 

equations , 


dN 


dt1  =  (n2  -  gNj)R21  +  n2c21  +  n5c51  +  Ne(N3K31  +  n4k41) 


+  N  3K*  -  N.  N  I  K.  +  K*  }] 

e  cl  1L  e 1  Ln  lm  lcJ  J 
m=2 


(4) 


dN, 


dt  =  ^gNl  "  N2^R21  +  Ne^NlK12  +  N5K52^  +  N3C32  +  N4C 
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*  -  N2[Ne{  I  K  *  K*c|  *  C21]  -  N^v(aLF  *  „A)  - 

m=3 


(5) 


and 


*In  order  to  obtain  the  close  agreement  (indicated  later)  between  the  5  and 
20-level  LIBORS  codes  it  was  found  necessary  to  set  Knc  =  0,  for 

n  =  17,  12  and  8,  i.e.,  the  6f,  5f  and  4f  levels  of  sodium.  As  a 

consequence,  p(n,  4)  =  1  for  n  =  17,  12  and  8  because  these  levels  are  only 

optically  coupled  to  level  4  of  the  manifold. 
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M  i  <NlKlm  *  N2K2n,  *  Ne2K™>)  -  N3<C32  +  NeK31  *  P3c> 
m=3 


-  N4(C42  +  NgK41  +  P4(;)  -  N5(C51  +  NgK52  +  P5c) 


where  we  have  introduced  the  (3-4-5)  level  manifold  population  density 


N*  =  N3  +  N4  +  N& 


(7) 


As  indicated  we  assume  that  the  close  spacing  of  these  three  levels  ensures 
collisional  equilibrium  and  consequently  that  the  internal  distribution  of 
this  manifold  is  described  by  a  Boltzmann  distribution.  In  which  case 


where 


N.  =  N*  /U . (T  ) 

J  3  j  '  e ' 

-E  ,/kT 

MV  =  l  9mie  mJ  6 

J  e  m=3,4,5  ^ 


(8) 

(9) 


with  g  .  =  g  /g .  and  E  .  -  E  -  E .  (j  =  3,  4  or  5) .  In  order  to  streamline 
3mj  3m  mj  m  j 

the  equations  as  far  as  possible  we  also  introduced  the  notation. 


and 


C  =  N  K  +  A 
nm  e  nm  nm 


P  =  a  F  +  N  K* 
nc  nc  e  nc 


(10) 
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where  a  is  the  single  photon  ionization  cross  section  for  level  n  and  F 
is  the  laser  photon  flux  density,  i.e.. 


F(t)  =  ^  /  I*(v,t)dv 


(12) 


where  I*(v,t)  is  the  laser  spectral  irradiance  at  time  t  and  hv  the  laser 
photon  energy.  In  the  above  equations  we  have  also  introduced  the 


resonance  to  ground  level  degeneracy  ratio  y  and  the  stimulated  emission 
rate  per  resonance  state  atom, 


g 

R2i (t )  =  i£r  I  IA(v,t)L21(v)dv  (13) 

where  is  the  Milne  stimulated  emission  coefficient  and  l_2^(v)  is  the 
line  profile  function  (we  shall  further  discuss  this  later)  for  the 
resonance  transition. 

The  appropriate  electron  energy  equation  can  be  expressed  in  the 
following  form: 


5  2  5  5 

4r(Ne)=N  V  Y  N  K  E  +  N  3  £  K*  E#  +  T  (E01  -  E  )N  o  F 

dt  '  e  e  e  n  nm  nm  e  L.  cm  cm  L~  v  21  cm7  m  me 
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3 

where  eg  represents  the  free  electron  mean  thermal  energy  (taken  as  ^-kTe), 
and  Hg. ,  Hgl  and  Hg2  represent  the  rates  of  elastic  energy  transfer  to  ions 
and  the  ground  and  resonance  state  atoms  through  Coulomb  scattering 
collisions^)  and  hard  sphere  collisions, respectively.  S  corresponds 
to  the  source  of  electron  energy  arising  from  the  various* seed  ionization 
processes,  i.e.. 


S  =  (2E21  '  Ec2)tN2G2c)p2  +  7  N22voLF)  + 
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Here  o2c  represents  the  resonance  state  two  photon  ionization  coefficient, 
the  laser  assisted  Penning  ionization  coefficient,  and  the  effective 
associative  ionization  cross  section  which  is  taken  to  be  an  F-dependent 
empirical  fit  to  the  experimental  work  of  Roussel  et  al.^^  E.  is  the 
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effective  energy  of  creation  for  the  electrons  liberated  by  this  latter 
process  and  v  represents  the  mean  atom  collisional  velocity.  Equation  (14) 
is  solved  in  conjunction  with  the  ionization  (Ng)  equatiion.  A  detailed 
discussion  of  the  20-level  atomic  code  and  a  comprehensive  review  of  the 
various  seed  ionization  processes  with  the  relevant  cross  sections  used  is 

( 3i ) 

provided  by  Cardinal.' 

It  should  be  noted  that  another  fine  adjustment  in  the  form  of  a  small 

reduction  (of  about  0.055  eV)  in  the  three-body  recombination  energy  was 

found  to  bring  the  results  of  the  5-level  LIBORS  code  into  excellent 

(251 

agreement  with  the  20-level  code  over  a  wide  range  of  densities.  '  This 

M 

modified  energy  is  designated  by  Ecm i n  equation  (14). 

The  electron  - ion  energy  exchange  rate  H  .  is  proportional  to 

(301  e 1 

(T  -T^)'  '  and  therefore  it  is  necessary  to  evaluate  the  ion  temperature. 

This  is  accomplished  through  simultaneously  solving  the  ion  energy 

equation , 


(fF  <Vl>  ■  Ne  Hei  -  (No  -  Wt.  *  WlA  *  <16> 
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where  we  assume  that  Ng  =  N.  and  and  Efl  are  the  mean  thermal  energies  of 

the  ion  and  atom,  respectively.  The  second  and  third  terms  in  equation 

(16)  takes  account  of  elastic  cooling  collisions  between  the  sodium  ions 

and  the  background  of  sodium  and  argon  buffer  gas  atoms.  (Nq  -  Ng )  and 

represent  the  sodium  and  argon  atom  densities  while  Hia  and  H^A  represent 

respective  rates  of  elastic  energy  exchange  between  the  ions  and  the 

(32  331 

background  of  neutral  sodium  and  argon  atoms'  *  ' 


where  <ov>  is  the  effective  ion-argon  collision  frequency  and  T  the 

1  A  A 

argon  atom's  translational  temperature.  A  similar  expression  is  used  in 

the  case  of  H.  ,  with  m.  =  m  in  this  instance.  The  last  two  terms  in  Eq. 
i  a  i  a 

(16)  takes  account  of  the  net  gain  or  loss  of  energy  to  the  ions  resulting 
from  their  creation  from  atoms  or  their  conversion  back  to  atoms. 
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This  background  density  of  argon  atoms  varies  across  the  sodium  vapor 

region  such  that  the  sum  of  the  partial  pressures  equals  the  pressure  of 

the  argon  in  the  water  cooled  peripheral  region.  The  initial  vapor 

temperature  at  each  location  is  deduced  from  the  optically  measured  sodium 

( 34) 

density  profile  using  the  vapor  pressure  curves  provided  by  Nesmeyanov.'  ' 

By  a  similar  argument  to  that  given  above  it  follows  that  since  H.^is 
proportional  to  the  difference  in  temperature  between  the  sodium  ions  and 
the  argon  atoms,  the  latter  shall  be  computed.  We  have  assumed  that  the 
argon's  temperature  is  not  changed  appreciably.  It  should  also  be  noted 
that  an  electron-argon  elastic  energy  loss  term  was  introduced  into  the 
free  electron  energy  equation.  To  complete  the  set  of  energy  equations  we 
have  also  included  an  equation  for  the  mean  thermal  energy  of  the  neutral 
sodium  atoms. 

Combining  equations  (12)  and  (13)  enables  us  to  write  the  stimulated 
rate  of  emission  per  atom  in  level  2  in  the  form 

2  oo 

R21(t)  =  rTT-  F(t)  f  L*(x  -x)L21(x)dx  (18) 

where  \  is  the  appropriate  wavelength  and  the  radiative  lifetime  of 
level  2  and  x  is  a  dummy  frequency  variable  introduced  below  with  x^ 
corresponding  to  the  laser  frequency. 

In  reality  level  2  for  the  sodium  atom  comprises  the  two  fine 
structure  levels:  32P3/2  and  32P1/2  and  so  the  reciprocal  of  the  effective 
1  i  fetime 

— -  =  j  [A(32P3/2,  32S1/2)  +  A(32P1/2,  3231/2)  ].  (19) 

T21  c 

Half  the  sum  of  the  Einstein  spontaneous  transition  probabilities  for  the 
two  resonance  transitions  is  taken  because  in  equations  (4)  and  (5)  we  use 
N^,  the  total  population  of  the  32P  term. 

The  corresponding  doublet  line  profile  function  for  the  sodium 
resonance  lines, 

^"2  1  ^  x  )  =  [‘3^3/2^x^  +  3  L 1  /  2(  X  +  A)]  (20) 


E.  10 


The  dummy  frequency  variable 


x  =  V  -  v3/2  (21) 

where  v3/2  represents  the  center  line  frequency  for  the  (32P3/2  -  32Si/2) 
transition  (at  589.0  nm)  and  A  =  v3/2  -  v1/2  is  the  frequency  separation  of 
the  doublet.  3 / x )  anc*  *-1/2(x+a)  are  the  respective  line  profiles  of  the 
individual  resonance  transitions.  We  have  assumed  a  Lorentzian  line 
profile  function  of  the  form 


L  (v  -  v  ).i - - -  (22) 

3/2  3/2  *  (v-v3/2)2  +  ya2 

where  is  the  appropriate  HMHW  bandwidth  [a  combination  of  the  Stark 
width  and  the  resonance  width], (25.26,35) 

We  have  also  introduced  the  spectral  distribution  function  for  the 
laser  radiation 


_  I  (v,  z,  t) 

‘  v)  ‘  "Wfz.  t) 


and  assumed  that  it  is  approximated  by  a  Gaussian,  i.e.. 


(23) 


Ll(vrv) 


-[(v.-v)/PJ2 
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(24) 


In  reality,  there  would  be  some  degree  of  spectral  hole  burning  expected  as 

the  laser  pulse  propagates  through  the  vapor  but  as  shown  by  Measures  and 
(36)  (26) 

Herchen'  ;  Wong,v  ;  the  omission  of  this  feature  does  not  lead  to 
appreciable  error.  Consequently,  the  convolution  of  the  two  profiles  in 
(18)  can  be  expressed  in  the  form  of  a  Voigt  function, 

2 

V(  6  ,u )  =— £ -  f  — -  (25) 
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71 

where 
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y  =  (x^-xj/p^. 


u  s 


and 


V(6,  u)du  =  V(6,  v)dv 


6  E  W  -  x  x  x  x 

It  is  important  to  notice  that  this  Voigt  function  depends  upon  6  -  the 
ratio  of  the  atomic  bandwidth  to  the  laser  bandwidth  and  to  the  frequency 
separation  of  the  atomic  and  laser  center  frequencies  normalized  by  the 
laser  bandwidth ,  i .e. , 


V*  "  V3/2 


The  form  of  equation  (18)  leads  us  to  introduce  an  "effective 
cross-section"  for  stimulated  emission  of  the  resonance  doublet, 

°21  =  — -  /  L*(x!  "  x)l_21(x)dx  (27) 

8ht21  -® 

We  recognized  from  the  outset  that  attenuation  of  the  laser  pulse 
could  lead  to  relatively  sudden  changes  in  the  laser  irradiance  and  that  in 
many  experiments  short  laser  pulses  are  used.  Consequently  in  our 
comparison  of  the  5  and  20-level  computer  codes  we  used  the  truncated  laser 
pulse  indicated  as  F  in  figure  4.  With  a  laser  pulse  of  this  shape  the 
comparison  was  tested  fairly  severely.  Figure  4  indicates  the  excellent 
agreement  obtained  between  the  5  and  20-level  models  for  a  sodium  density 
(N0)  of  4*1016  cm-3.  It  is  evident  from  figure  4  .that  N2/Nx  falls  from  the 
saturated  value  of  3.0  at  about  the  time  of  ionization  burnout  even  though 
F  is  still  increasing.  This  temporal  minimum  is  attributed  to  the  sudden 
increase  in  the  superelastic  quenching  rate  which  for  the  given  laser 

irradiance  can  overwhelm  the  stimulated  emission  rate  as  full  ionization  is 
achieved. 

It  is  apparent  from  figure  4  that  the  5-level  LIBORS  code  mirrors  very 
closely  the  predictions  of  the  20-level  code  in  regard  to  the  temporal 
behaviour  of  the  parameters:  Ne>  Tg,  T^ ,  and  N^.  We  feel  that  the 

small  (<  20%)  difference  in  the  lead  time  between  the  two  sets  of  curves  is 

quite  acceptable.  In  order  to  check  that  this  close  agreement  was  not 
fortuitous  to  this  particular  density  we  ran  comparisons  at  two  additional 
densities  (Nq  =  5*1015  and  1016  cm-3)  and  found  the  same  close 


ayreenient . 
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RADIATIVE  TRANSFER  EQUATION  FOR  AN  EXTENDED  VAPOR  RE G I ON 


The  large  cross  section  involved  in  laser  resonance  saturation  leads 
to  appreciable  absorption  of  the  laser  pulse  as  it  propagates  through  the 
atomic  vapor.  Consequently,  any  realistic  simulation  of  plasma  channel 
formation  based  on  this  interaction  must  take  account  of  this  fact  and  we 
are  therefore  led  to  consider  the  radiative  transfer  equation.  For  a  laser 
pulse  propagating  in  the  z-direction  we  can  write 
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The  first  term  on  the  right  hand  side  takes  account  of  the  net  loss  of 
laser  photons  due  to  absorption  in  the  resonance  transition,  the  second 
term  allows  for  two  photon  ionization  of  the  resonance  level,  the  third 
term  involves  laser  assisted  Penning  ionization  of  the  resonance  level  and 
the  last  term  refers  to  single  photon  ionization  of  the  levels  within  the 
(3-4-5)  level  manifold. 

For  the  time  scale  of  interest  in  our  experiments  the  velocity  of 
light  c  can  be  treated  as  infinite.  Under  these  circumstances  (28) 
is  simplified  to  the  steady  state  form  given  by 


5 

f  *  <N2  -  -  N?42c  p2  -  7  -A  Vn,cF  <29> 


m=3 
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where  o21  was  given  by  equation  (27). 

It  is  quite  clear  that  addressing  the  problem  of  laser  resonance 
saturation  for  an  extended  vapor  region  involves  the  solution  of  the 
population  and  energy  equations  indicated  earlier,  together  with  this 
radiative  transfer  equation  in  both  space  and  time.  We  decided  to  tackle 
this  problem  by  dividing  the  irradiated  vapor  column  into  a  series  of 
slabs,  see  figure  5,  each  of  which  was  assumed  to  have  a  uniform  atom 
density  equal  to  the  average  experimentally  measured  sodium  atom  density 


for  that  segment  of  vapor.  The  thickness  of  each  slab  was  determined  by 

the  requirement  that  at  least  99%  of  the  energy  incident  on  any  yiven  slab 

(331 

be  transmitted  to  the  next  slabv  '  that  the  maximum  value  be  limited  to 
1  mm 

Within  each  slab  the  population  and  energy  equations  are  numerically 

(331 

integrated  (using  an  IMSL  subroutine  DGEAR'  '),  subject  to  the  laser  pulse 

exiting  the  previous  slab  to  provide  the  temporal  history  of  the  various 

parameters:  N  ,  T  ,  T.  Ni,  N0  and  N. .  In  this  multislab  model,  equation 
e  e  i  ^  i 

(29)  is  solved  and  the  laser  photon  flux  density  exiting  the  i  slab 
(bounded  by  z ^  and  z^ _j) 

F(z.,  t)  =  F(zi_1,  t)e  PiAZ  (30) 

and  is  seen  to  be  expressed  in  terms  of  the  incident  value  F(z^_^,  t).  In 
this  equation 


p.  =  [(gNj  -  N2)o21  + 


1 


>  & 

NoOo^F(z-  t)  +  4  N/vo.  +  l  N  a  1. 

2  2c  i-l  2  2  L  m  me Ji  (31) 
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and 

AZ  =  z.  -  (32) 

This  is  permissible  because  the  loss  of  laser  energy  through  two  photon 
ionization  is  small  compared  to  that  extracted  through  resonance  pumping. 
The  right-hand  side  of  equation  (31)  is  evaluated  using  the  mean  values  of 
the  parameters  associated  with  the  ith  slab.  The  value  of  F(z^,  t)  for 
each  time  step  is  stored  and  then  used  for  calculating  the  temporal 
variation  of  the  parameters  in  the  subsequent  (i+l)th  slab. 

The  incident  temporal  laser  profile  used  in  these  calculations  was 
chosen  to  approximate  the  experimental  laser  pulse  by  an  analytical  curve 
of  the  form 


F*(o,  t) 


.  b  -ct 
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(33) 


where  b  and  c  were  constants  determined  from  the  experiment  and  is  the 
eneryy  fluence  in  the  laser  pulse  ‘a"  is  determined  from  the  relation, 


r  *  b  -Ct2  ..  , 

J  at  e  dt  =  1 


EXPERIMENTAL  FACILITY 

The  experimental  facility  developed  for  this  research  is  schematically 
illustrated  in  figure  6.  The  second  harmonic  of  a  Neodymium-YAG  laser 
(JK  HY750)  was  used  to  pump  a  dye  laser  (Quanta-Ray  PD1).  Using  Kiton 
Red-620  dye  (Exciton  Chemical  Co.)  an  output  of  about  40  mJ  at  a  10  Hz 

repetition  rate  was  obtained  over  a  wavelength  range  that  included  the  two 
sodium  D  lines  at  589.0  and  589.6  nm.  The  dye  laser's  spectral  full  width 
half  maximum  FWHM  was  measured  to  be  about  0.02  nm. 

A  specially  designed  heat  sandwich  oven  was  used  to  generate  a 
disc-shaped  zone  of  sodium  vapor.  This  oven  is  in  the  form  of  a 

rectangular  box  with  circular  heater  plates  and  steel  wicks,  top  and 

bottom,  surrounded  by  cooling  water  coils  to  prevent  the  sodium  vapor  from 
reaching  the  four  windows  that  provide  360°  optical  access.  A  more 
detailed  description  of  the  method  of  measuring  the  sodium  atom  density 

/  nc  \ 

distribution  is  provided  elsewhere.  ' 

The  electron  density  measurements  were  undertaken  using  a  sideways 
mounted  SPEX  1700  II  monochromator.  This  enables  the  entrance  slit  of  the 
monochromator  to  sample  a  thin  horizontal  slab  of  the  plaspia  emission  at  a 
height  y  (scanned  by  adjusting  the  input  optics)  above  the  axis  of  the 
laser  beam.  The  magnification  of  the  optics  was  1:1  with  the  SPEX  entrance 
si  it  at  1  cm  x  40  pm. 

The  signal  from  the  RCA  C31034  photomultiplier  mounted  on  the  SPEX 
monochromator  was  processed  by  an  EG&G  4420  Signal  Averager  that  was 
typically  gated  to  sample  the  emission  in  a  2  ns  interval,  some  65  ns  after 
the  start  of  the  laser  pulse.  Since  the  laser  pulse  duration  was  less  than 
40  ns,  the  influence  of  the  laser  field  (through  the  dynamic  Stark  effect) 
should  be  minimal  and  spontaneous  and  electron  collision  induced  decay  of 
the  resonance  state  population  should  minimize  optical  depth  effects. 


In  order  to  evaluate  the  radial  profile  of  the  free  electron  density 
N  (4r),  spectral  scans  of  the  42D-32P  multiplet  were  taken  for  15  lateral 
(y)  positions  across  the  plasma  column.  The  spectral  resolution  was 
estimated  to  be  0.05  nm  (FWHW)  with  16  laser  shots  being  averaged  to 
constitute  one  intensity  measurement  (at  a  given  wavelength)  and  128  such 
measurements  comprise  one  spectral  scan. 

COMPARISON  OF  COMPUTATIONAL  AND  EXPERIMENTAL  RESULTS 

The  present  form  of  our  LIBORS  computer  code  predicts  the  state  of  the 
plasma  formed  along  the  direction  of  propagation  for  a  given  sodium  atom 
distribution  and  a  temporally  prescribed  laser  field.  In  an  experiment  the 
laser  field  also  has  a  radial  variation  and  this  has  to  be  taken  into 

account  if  we  wish  to  model  the  three  dimensional  nature  of  the 
interaction.  To  accomplish  this  we  have  assumed  that  the  incident  laser 
pulse  has  a  Gaussian  radial  distribution  (the  justification  for  this  is 

presented  el  sewhere ) ,  and  the  same  temporal  history  given  by  equation 
(33),  at  each  radial  position. 

A  series  of  computer  runs  for  a  range  of  incident  laser  energy 
fluences  was  then  undertaken.  The  state  of  the  plasma  at  any  position  can 
then  be  predicted  from  these  computer  runs  by  assigning  tne  appropriate 
laser  energy  fluence  to  each  radial  position.  By  way  of  example  we  have 

displayed  in  figure  7  the  predicted  axial  variation  in  the  free  electron 

density  Ng(z)  corresponding  to  four  radial  positions  (r  =  0»  1.25,  2.55  and 
3.55  mm)  65  ns  after  the  start  of  the  laser  pulse  for  the  experimentally 
based  sodium  atom  distribution  Nq(z)  also  shown.  The  incident  laser  pulse 
used  to  generate  these  results  wasa  assumed  to  have  the  temporal  variation 
given  by  equation  (33)  with  b  =  2.64,  c  *  6.5><10-3  and  =  40  ns  and  a 
Gaussian  radial  profile  with  an  energy  of  25  mJ  and  rQ  =  2.5  mm.  These 
four  radial  positions  corresponded  to  incident  laser  energy  fluences  of: 
127.3,  95,  55  and  17  (mJ  cm-2)  for  this  particular  laser  pulse. 

A  comparison  of  the  free  electron  density  achieved  at  equidensity 
axial  locations  in  figure  7  reveals  that  both  the  electron  density  and  the 
radius  of  the  plasma  column  co^e  diminishes  as  the  laser  pulse  penetrates 
further  into  the  sodium  vapor.  This  effect  is  even  more  obvious  in  figure 


8  where  the  predicted  full  radial  profile  of  the  electron  density  creott: 
by  this  laser  pulse  is  displayed. 

This  reduction  and  narrowing  of  the  core  plasma  density  is  a  direct 
consequence  of  the  absorption  suffered  by  the  laser  pulse  as  it  propagates 
through  the  sodum  vapor  (from  -ve  to  +ve  z  values).  In  figure  9  we  present 
the  axial  variation  of  the  fraction  of  transmitted  laser  energy  fluence  for 
several  incident  values  of  laser  energy  fluence  for  the  sodium  atom 
distribution  shown  in  figure  7.  It  can  be  seen  that  the  smaller  the 
incident  laser  energy  fluence  the  greater  its  percentage  attenuation  in 
propagating  through  the  sodium  vapor.  It  follows  that  the  radial  profile 
of  a  laser  pulse  will  tend  to  steepen  as  it  propagates  since  its  high 
intensity  core  will  be  proportionately  less  reduced  than  its  weaker  outer 
region . 

In  contrast  to  the  rapid  decline  of  t*v  core  electron  density  with 
axial  position,  the  low  electron  density  halo,  see  figure  8,  is  predicted 
to  fall  less  rapidly  with  axial  position.  This  behaviour  can  be  understood 
in  terms  of  the  shoulder  in  the  dependence  of  Ng  on  the  incident  laser 
energy  fluence  for  a  given  sodium  atom  distribution.  This  is  illustrated 
for  the  simulation  under  consideration  in  figure  10,  where  the  variation  of 
N  at  t  =  65  ns  and  Z  =  -4,  0  and  4  cm  are  plotted  as  a  function  of  the 
incident  laser  energy  fluence. 

We  can  interpret  the  N^-shoulder  seen  on  the  Z  -  4  an  curve  [30  c  < 

110  (mJ  cm-2)]  in  terms  of  our  four  stage  model  of  LIBORS  discussed 
(19) 

earlier.  '  For  >  110  (mJ  cm-2)  full  ionization  burnout  is  approached, 
while  in  the  range:  30  <  110  (mJ  cm-2)  the  degree  of  ionization 

attained  is  determined  by  the  energy  invested  in  saturating  the  atoms. 
This  resonance  state  energy  is  then  converted  into  ionization  through 
superelastic  collisional  heating  of  those  electrons  produced  by  the  seed 
ionization  process.  In  essence  each  atom  can  only  contribute  one  laser 
photon  in  this  regime  (compared  to  being  cycled  many  times  with  longer  more 
intense  laser  radiation)  and  so  the  degree  of  ionization  attained  is  almost 
independent  of  the  laser  energy  fluence  under  these  conditions.  This 
suggestion  is  borne  out  when  the  computer  predicted  energy  loss  per  cm  path 
length  is  compared  with  the  energy  required  to  saturate  the  atoms  per  cm  at 
this  location. 
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In  figure  11,  we  display  the  experimentally  measured  free  electron 

density  radial  profiles  at  Z  =  -2  and  2  cm,  65  ns  from  the  start  of  the 

incident  laser  pulse.  These  results  were  derived  by  fitting  convoluted 

Stark  and  instrumental  profiles  to  radially  inverted  using  the  Abel 

transform  42D-32P  lateral  emission  spectra^)  for  the  same  conditions  as 

used  in  the  computational  work.  Also  shown  in  figure  11  are  two  empirical 

profiles  fitted  to  these  data.  A  comparison  of  figure  11  with  that  of 

figure  9  reveals  that  the  experimental  Ng-radial  profiles  fall  off  more 

sharply  than  the  predicted  radial  profiles.  In  light  of  the  above 

discussion  this  led  us  to  rerun  our  computer  code  with  a  reduced  value  for 

the  associative  ionization  cross-section  o^.  To  simplify  matters  we  used  a 

constant  value  a.  =  1.5xl0-15  cm2  rather  than  scale  the  irradiance 

”  fill 

dependent  curve  suggested  by  Roussel  et  al.'  '  This  is  justified  for  our 
computations  since  had  its  maximum  value  of  5*10“15  cm2  for  most  of  the 
conditions  of  interest,  the  exception  being  large  radial  positions  (r  >  3 
mm)  at  axial  locations,  Z  >  2  cm. 

In  figure  12(a)  we  present  both  the  experimental  N  -radial  profile  (at 

C 

Z  =  -2  cm  and  t  =  65  ns )  and  the  two  computed  profiles:  for  o.  =  o.(F), 

/  i  i  \  MM 

derived  from  Roussel  et  al,  '  and  oA  =  1.5x10” 15  cm2,  while  in  figure 
12(b)  we  present  the  corresponding  profiles  at  Z  *  2  cm.  It  is  evident 
that  the  experimental  values  of  Ng(r)  fall,  to  some  extent,  between  the  two 
computational  curves  at  both  Z  =  -2  and  2  cm.  The  experimentally  derived 
radial  profile  is  seen  to  be  steeper  than  either  of  the  computational 
profiles,  in  keeping  with  our  earlier  observation.  Indeed  both  computed 
profiles  underestimate  the  peak  ionization  and  overestimate  the  ionization 
at  large  radial  positions.  We  should,  however,  point  out  that  our  sodium 
atom  density  measurements  are  only  accurate  to  about  20%^^^  and  that  our 
assumed  radius  of  the  laser  pulse  also  is  subject  to  comparable 
uncertainty.  We  will  address  this  further  in  the  Discussion  and  Conclusion 


section. 

In  addition  to  the  radial  measurements  of  the  free  electron  density  we 
have  also  recorded  the  42D-32P  multiplet  spectrum  at  five  (y=0)  axial 
locations  (Z  =  -4,  -2,  0,  2,  4  cm).  These  spectra  can  be  used  to  give  us  a 
fairly  good  idea  of  the  variation  of  the  peak  (r  =  0)  electron  density 
along  the  path  of  the  laser  beam.  This  is  accomplished  by  matching  the 
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spectrum  at  each  Z-location  to  that  evaluated  from  the  one  dimensional 

radiative  transfer  equation  for  a  uniform  plasma  of  electron  density  N  and 

(371  e 

radius  rQ.v  '  In  effect  Ngand  rQ  are  used  as  fitting  parameters.  Figure 

13  presents  these  two  sets  of  42D-32P  multiplet  spectra.  The  computed 

Stark  broadened  profiles  (full  curves)  are  seen  to  approximate  the 

experimental  spectra  (+)  in  regard  to:  the  ratio  of  peak  heights,  the  peak 

to  minimum  ratio  and  the  widths  of  the  profiles. 

Confirmation  that  these  spectra,  which  are  integrated  along  the  line 
of  sight  for  y  =  0  (as  opposed  to  the  radially  inverted  spectra  used  for 
the  radial  profile  measurements)  can  be  used  to  approximately  assess  the 
axial  (r  =  0)  electron  density  is  obtained  when  we  compare  the  predictions 
at  Z  *  -2  and  2  cm  from  figure  13  with  the  actual  radial  profiles  presented 
as  figure  11. 

The  axial  values  of  Ng  can  be  seen  to  decline  rapidly  for  Z  >  0  and 
the  asynroetry  with  regard  to  equidensity  locations  is  in  keeping  with  our 
other  observations  and  with  our  computational  predictions.  These  results 
confirm  that  attenuation  of  the  laser  pulse  is  a  very  real  effect,  even  for 
distances  of  a  few  centimeters,  once  the  atom  density  is  around  10 16  cm-3. 

If  we  assume  that  the  plasma  is  in  LTE,  then  the  experimentally 
measured  radial  electron  density  distributions  at  Z  *  -2  and  2  cm  can  be 
used  with  the  Saha  equation  and  the  initial  atom  density  to  determine  the 
corresponding  radial  profile  of  the  free  electron  temperature.  These  are 
displayed  with  the  two  computed  Tg(r)  profiles  in  figure  14(a)  and  (b) .  In 
both  cases  the  experimental  data  is  quite  well  flanged  by  the  two 
computational  T  -profiles. 

The  temperatures  indicated  in  these  results  are  somewhat  lower  than  we 

had-  previously  expected  from  LIBORS  but  appears  to  confirm  other 
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observations  of  low  electron  temperatures.  ’  ’  This  is  particularly 

so  for  Z  =  2  cm,  where  the  mean  temperature  is  about  4000K.  These  low 
temperatures  are  predicted  by  our  LIBORS  code  after  penetrations  of  several 
centimeters  of  vapor  and  represent  another  manifestation  of  reduced  laser 
energy  fluence  resulting  from  absorption  of  the  laser  pulse.  This  is 
clearly  seen  in  figure  15  where  the  computer  predicted  axial  variation  of 
Ng  and  Tfi  (for  r  *  0)  are  presented  with  that  of  for  =  1.5x10" 15  cm2. 
Also  shown  is  the  measured  axial  variation  of  NQ  that  was  used  in  the 
computer  code.  The  slight  upturn  in  Tg  predicted  at  Z  *  -5  and  4  cm  is  a 


result  of  the  ionization  not  proceeding  to  the  burnout  stage.  The 
reliability  of  the  computer  calculations  is  questionable  at  these  low 
densities  since  the  free  electrons  may  not  Maxwellianize  in  the  time 
allowed.  In  figure  16  we  present  a  comparison  of  the  computer  predicted 
axial  variation  (at  r  =  0)  of  Ng  and  Tg  with  that  of  the  experimentally 
derived  values. 

The  cmputational  work  used  the  associative  ionization  cross  section  of 
Roussel  et  al^1^  and  is  found  to  give  results  that  are  in  close  agreement 
with  the  experimental  values.  The  deviation  observed  for  Z  >  0  could  be 
due  to  the  uncertainty  in  the  sodium  atom  density  distribution.  A  slightly 
large  Nq (z )  would  increase  the  computed  Ng  values  for  Z  <  0,  but  lower  them 
for  Z  >  0  due  to  attenuation  of  the  laser  pulse. 


DISCUSSION  AND  CONCLUSIONS 

We  have  developed  a  computational  model  of  "laser  ionization  based  on 
resonance  saturation"  (LIBORS )  for  sodium  vapor  that  takes  into  account 
both  the  nonuniform  atom  distribution  and  absorption  of  the  laser  field. 
This  computer  code  has  been  used  to  provide  the  first  insight  into  the 
three  dimensional  nature  of  this  interaction. 

We  have  also  undertaken  the  first  measurements  of  both  the  radial  and 
axial  profiles  of  the  free  electron  density  in  a  sodium  plasma  created  by 
laser  resonance  saturation.  From  these  measurements  we  have  deduced  the 
corresponding  radial  and  axial  electron  temperature  profiles. 

We  have  used  our  LIBORS  cmputer  code  to  model  the  sodium  plasma  formed 
along  the  path  of  the  laser  pulse  used  in  our  experiments  and  have  compared 
the  code  predicted  radial  and  axial  profiles  of  electron  density  and 
temperature  with  the  results  derived  from  our  experiments.  In  general  the 
agreement  is  very  good,  especially  if  allowance  is  made  for  the 
experimental  uncertainties  and  the  limited  accuracy  of  the  many 
cross-sections  used  in  the  code. 

A  detailed  comparison  reveals  that  the  computer  code  (using  o.  of 

(111  ” 
Roussel  et  al v  '  tends  to  underpredict  the  r  =  0  value  of  the  free 

electron  density  for  Z  <  0,  while  overpredicting  the  value  for  Z  >  0.  This 

is  consistent  with  a  few  percent  underestimate  of  the  sodium  atom  density. 

This  would  lead  to  larger  values  of  Ng  for  Z  <  0,  but  lower  values  of  Ng 

for  Z  >  0  due  to  greater  attenuation  of  the  laser  pulse.  The 


mmmm 
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overprediction  of  the  code  with  respect  to  the  free  electron  density  at 
large  radial  positions  could  be  accounted  for  by  an  overestimate  on  our 
behalf  of  the  radius  of  the  incident  laser  pulse. 

We  believe  that  these  results  demonstrate  that  our  LIBORS  computer 
code  is  capable  of  fairly  accurately  predicting  the  three  dimensional 
behaviour  of  this  new  mode  of  laser  ionization.  In  particular,  both  our 
computer  code  and  our  experiments  clearly  reveal  that  absorption  of  laser 
energy  due  to  the  strong  nature  of  the  interaction  leads  to  a  decline  in 
the  free  electron  density  and  temperature  along  the  path  of  the  laser 
pulse.  This  can,  for  modest  values  of  laser  energy  fluence,  lead  to  quite 
low  values  of  both  Ng  and  T0  after  the  laser  pulse  has  penetrated  just  a 
few  centimeters  of  vapor  at  an  atom  density  of  about  1016  cm-3.  This  could 
well  explain  why  some  researchers  have  failed  to  observe  appreciable 
ionization or  have  measured  low  values  of 

Our  results  clearly  indicate  that  if  LIBORS  is  to  be  used  to  generate 

fairly  uniform  plasma  channels,  then  the  laser  energy  fluence  has  to  be 

sufficient  to  achieve  close  to  full  ionization  along  the  path  of  the  atomic 

(41  42) 

vapor.  If  plasma  channels  of  several  meters  length  are  required,'  * 
this  could  require  very  large  pulses  of  laser  energy,  and  other  techniques 
such  as  near  resonant  multiphoton  ionization  could  be  quite  competitive, 
even  if  they  do  not  achieve  full  ionization  during  the  laser  pulse 
duration. 

In  comparing  LIBORS  with  near  resonant  multiphoton  ionization  it  has 
to  be  recognized  that  in  the  former  case  ionization  is  achieved  through  a 
laser  driven  collisional  process  which  invariably  leads  to  populating  most 
of  the  intermediate  levels  —  all  of  which  lose  energy  by  radiative  decay. 
Near  resonant  multiphoton  ionization  can  lead  to  direct  ionization  with 
almost  no  excitation  of  intermediate  states  for  short  times.  However,  it 
is  difficult  to  achieve  full  ionization  by  this  means  especially  at  high 
densities.  Consequently  LIBORS  may  be  the  best  technique  for  creating 
long,  dense,  relatively  cool  plasmas. 
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Fig.  16  Comparison  of  the  axial  variation  (for  r  =  0)  of  the  free  electron 
density  and  temperature  determined  by  experiment  and  predicted  by 
our  LIBORS  computer  code. 
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APPENDIX  F 


SUPERELASTIC  COLLISIONAL  EXCITATION  OF  A  STRONTIUM  PLASMA 
BY  ONE  AND  TWO  PHOTON  LASER  PUMPING 

K.  A.  Menard  and  R.  M.  Measures 

Institute  for  Aerospace  Studies,  University  of  Toronto 
4925  Dufferin  Street 
Downsview,  Ontario,  Canada 
M3H  5T6 


Abstract 

Fast  superelastic  electron  heating  by  collision  quenching  of  one  and 
two  photon  laser  pumped  states  has  been  demonstrated  in  a  relatively  cold 
strontium  plasma  created  by  two  photon  ionization  of  strontium  vapor 
contained  within  a  crossed  heat  pipe  oven. 


nit,  under  certain  circumstances,  this  form  of  laser  Interaction  could  be  and  a  s ) i ynt 1 y  stronger  signal  at  the  b2Fi/2,//2  *  s,2bs/2  wavelength 
iery  competitive  with  inverse  bremsstrahl ung  for  coupling  laser  energy  into  (B6B.2U  to).  To  study  the  sensitivity  of  this  excitation  process  to  laser 
i  plasma.  Measures  et  al^21^  also  suggested  that  multiphoton  pumping  could  wavelength,  the  signal  was  temporally  averaged  over  the  length  of  the  laser 
>e  used  to  broaden  the  scope  of  this  form  of  laser  interaction.  We  suspect  pulse,  using  the  gated  sampling  mode  of  the  Boxcar  Averager  and  the  laser 


Figure  1.  Schematic  View  of  Strontium  Multiphoton  Pumping  Experimental 
Facility. 
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four-wave  difference  mixing  has  been  demonstrated  in  any  ion  using  a  single 
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oemberyen,  this  is  a  very  Interesting  result.  Even 


,he  phot odetect ion  system  comprised  an  RCA  4526  photomultiplier  and  transition,  the  luminous  intensity  Increased  and  shifted  to  a 

•'tral  discrimination  was  accomplished  with  a  SPEX  1700-11  Cterny-Turner  predominately  blue-white  colour.  The  most  intense  breakdown  occurred 


